m 

o 

(N 

< 

(N 






I 

o 
o 



m 

cn 

o 
en 



% 



Space group symmetry, spin-orbit coupling and the low energy effective Hamiltonian 

for iron based superconductors 

Vladimir Cvetkovic^ and Oskar Vafek^ 

^National High Magnetic Field Laboratory and Department of Physics, 

Florida State University, Tallahasse, Florida 32306, USA 

(Dated: April 16, 2013) 

We construct the symmetry adapted low energy effective Hamiltonian for the electronic states 
in the vicinity of the Fermi level in iron based superconductors. We use Luttinger's method of 
invariants, expanding about V and M points in the Brillouin zone corresponding to two iron unit 
cell, and then matching the coefficients of the expansion to the 5- and 8-band models. We then 
use the method of invariants to study the effects of the spin-density wave order parameters on 
the electronic spectrum, with and without spin-orbit coupling included. Among the results of this 
analysis is the finding that the nodal spin-density wave is unstable once spin-orbit coupling is 
included. Similar analysis is performed for the Aig spin singlet superconducting state. Without 
spin-orbit coupling there is one pairing invariant near the F point, but two near the M point. This 
leads to an isotropic spectral gap at the hole Fermi surface near F, but anisotropic near M. The 
relative values of these three parameters determine whether the superconducting state is S4.+ , s+_, 
or nodal. Inclusion of spin-orbit coupling leads to additional mixing of spin triplet pairing, with 
one additional pairing parameter near F and one near M. This leads to an anisotropic spectral gap 
near both hole and electron Fermi surfaces, the latter no longer cross, but rather split. 



I. INTRODUCTION 

The parent state of most iron based superconductors is 
an itinerant spin-density wave which, upon dopingir— or 
pressure^, gives rise to superconductivity. Unlike their 
copper oxide counterparts, which are doped Mott insula- 
tors with an odd number of electrons per unit cell, such 
iron based superconductors arc compensated semi-metals 
with an even number of electrons per unit celliSiii. The 
correlation effects in the iron based superconductors ap- 
pear to be more significant than in conventional metals, 
but not as severe as in copper oxide superconductors^^. 
Unlike in conventional superconductors, the mechanism 
of superconductivity is believed to be different from the 
electron-phonon interaction driven pairingi^^— . 

This paper is motivated by the need to develop a low 
energy effective theory which successfully describes the 
electronic states in the vicinity of the Fermi level. Such 
low lying states which are responsible for a number of 
physical characteristics of these materials, such as elec- 
trical and thermal conductivity, low temperature specific 
heat and magnetic susceptibility, are most affected by 
spin or charge ordering, or by superconductivity. Cur- 
rently, the models which successfully capture these states 
are either based on the 5-band tight-binding approach 
or by starting with it, diagonalizing, and then working 
in the resulting "band basis" as if it were a continuum 
model. The former has the disadvantage of being im- 
practical for studying the effects of an externally ap- 
plied magnetic field, or, simply, because it necessitates 
working with a large matrix. The latter also cannot be 
minimally coupled to the external magnetic field. More- 
over, it has an additional disadvantage, which stems from 
the k-space degeneracy of the hole bands at the F-point 
and the electron bands at the M-point, and which re- 
sults in k-space non-analyticity of the single particle 



wave-functions. Such non-analyticity translates into non- 
local effective electron-electron interaction, making the 
method impractical for studying interaction effects. 

In contrast, the approach which we develop here al- 
lows minimal coupling to the external magnetic vector 
potential, and maintains locality of the electron-electron 
interactions, provided they start out being of finite range 
in the lattice model. The method is based on the the- 
ory of invariants used by Luttingei— to study Si and 
Ge. The connection to the methods used to study 
semiconductorslSi is natural given that the parent state of 
iron-based superconductors is a multi-band semi-metal. 
The interesting new aspect is the manifest presence of a 
broken symmetry, which can be readily included within 
this approach. Furthermore, the space group of the iron- 
pnictogen or iron-chalcogen plane (we use "iron plane" 
from now on) is non-symmorphic because it contains an 
?i-glide plane: a mirror reflection about the iron plane 
followed by the translation along the half of the unit cell 
diagonal. This has profound consequences on the na- 
ture of the irreducible representations of the group of 
the wave- vectors at the Brillouin zone edges or corners, 
precisely where the electron pockets appear— i^Sr—. As 
we show below, the group of the wave-vector at M has 
only two dimensional irreducible representations, which 
means that, not including spin, at M all Bloch states 
are doubly degenerate. This is unlike k-points inside the 
Brillouin zone, such as F-point where the hole pockets 
appear, where the irreducible representations correspond 
to a known point group and therefore always have one di- 
mensional representations — if such point group is non- 
abelian, then of course it also has higher dimensional, ir- 
reducible representations, but it always has some one di- 
mensional irreducible representations. The states which 
cross the Fermi level near F derive from the Fe dxz and 
dyz orbitals and transform under the Eg representation of 



the group of the wave- vector which is isomorphic to D4/1. 
This representation is two dimensional, guaranteeing the 
degeneracy of these two states above any structural or 
magnetic transition. 

The method used here allows us to analyze the effects 
of the atomic spin-orbit interaction, AL-S, which, in iron, 
has been rcported^S to be A ~ 8QmeV. Such a value is 
larger than the temperature scale associated with both 
the magnetic and superconducting ordering. We find 
that once the spin-orbit coupling is included, one cannot 
avoid using the two iron unit cell: the two electron Fermi 
surfaces near the M-point no longer intersect, but rather 
split due to avoided level crossing. Moreover, the degen- 
eracy at the F-point is lifted through a spin-orbit term 
which is analogous to the one written by Kane and Mele 
in graphene^. When the time reversal symmetry and 
a center of inversion are present, as in the normal non- 
magnetic state, all bands are Kramers degenerate. At 
the M-point, the states are four- fold degenerate, due to 
a single four dimension double-valued irreducible repre- 
sentation of the space group*^. In the immediate vicinity 
of such four-fold degeneracy, the bands disperse linearly 
in momentum, as for a massless Dirac particle. For larger 
momentum deviation, the bands disperse upward. 

We also analyze the electronic spectrum in the pres- 
ence of the coUinear spin density wave^i^. Such a state 
breaks time reversal and inversion symmetry, but not 
their product. Therefore, the bands remain Kramers de- 
generate even below the spin-density wave ordering tran- 
sition temperature. For weak spin-density wave order, 
we find that the electron and hole Fermi surfaces recon- 
struct, leaving behind several smaller pockets. In the 
absence of spin-orbit coupling, we find 6 nodal points be- 
low the Fermi level, similar to Y. Ran et.al^. However, 
once the spin-orbit interaction is included, the degener- 
acy at the nodal points is lifted. This can be understood 
as a consequence of (a generalization of) the Wigner- 
von Neuman argument for Kramers degenerate bands. 
As a result, we conclude that for a putative strong spin- 
density wave, where Y. Ran et.alM predicted the gapless 
spectrum with nodal (Dirac) points, the spin-orbit cou- 
pling results in a fully gapped spectrum. In the case of 
co-planar spin-density wave with four-fold symmetry, re- 
cently reported in Ref. 1351 . the center of inversion is lost 
and the Kramers degeneracy is lifted. The resulting spec- 
trum therefore displays split Fermi surfaces. Based on 
the symmetry invariants that can be constructed, one can 
readily see that a transition from a normal non-magnetic 
state into a collincar magnetic state necessarily induces 
orthorhombic distortion, but not vice versa, in agree- 
ment with, for example, Refs. 13614401 A transition from 
a coUinear spin-density wave into a coplanar spin-density 
wave with 4-fold symmetry is necessarily of first order, in 
agreement with Ref. [3^. The absence of the Kramers de- 
generacy in the coplanar spin-density wave state suggests 
interesting repercussions on the microscopic coexistence 
of such a state with superconductivity, similar to Ref. l4ll . 

This method also allows us to analyze the symmetry of 



the superconducting pairing states and the Fermion bilin- 
ear terms they induce in the Hamiltonian. In the absence 
of spin-orbit coupling, at the F point there is a single, mo- 
mentum independent, spin singlet Aig invariant, whose 
strength we parameterize by Ar, which apart from charge 
L'^(l), does not break any other symmetry of the crystal. 
This term leads to an isotropic spectral gap along the 
two (anisotropic) hole Fermi surfaces. Similarly, at the 
M point, there are two such invariants, whose strength 
we parameterize by Ami and Ams- Depending on the 
relative ratio of Ami and Am3 we find the spectrum to 
be either gapped or nodal at the electron Fermi surfaces. 
Once spin-orbit coupling is included, there is an addi- 
tional, momentum independent, invariant at the F point. 
We parameterize it by Art. While being spin triplet, 
it nevertheless respects all the symmetries respected by 
Ar. The inclusion of this term makes the spectral gap 
on the hole Fermi surfaces anisotropic. At the M-point, 
the spin-orbit coupling also permits one additional mo- 
mentum independent triplet invariant, which contributes 
to the anisotropy of the spectral gap along the two, now 
split, Fermi surfaces. 

This paper is organized as follows: in the next subsec- 
tion we summarize the key results and point the reader 
to the equations corresponding to them. In Section [TTl 
we present an extensive analysis of the irreducible repre- 
sentation of the space group of iron based superconduc- 
tors, provide the product tables of the irreducible rep- 
resentations, illustrate the construction of the symmetry 
adapted functions at the M-point, and analyze the sym- 
metry of the Bloch states in the vicinity of the Fermi 
level. In Section III, we construct our low energy effec- 
tive Hamiltonian in terms of the six component Fermi 
"spinor" and determine the coefficients of the symmetry 
allowed terms to the 5-band tight-binding models of Ref. 
I43 and the 8-band model of Ref. [III. The resulting Hamil- 
tonian is very easy to handle and effectively requires di- 
agonalization of only a 2 x 2 matrix. We also include the 
effects of spin-orbit coupling and critically compare our 
approach to the two- and three-orbital models used in the 
literature. In addition, we use our machinery to build 
the symmetry allowed four-Fermion contact interaction 
terms. In Section HVl we study the consequences of the 
spin-density wave ordering on the electronic spectra, and 
similarly in Section V we study the superconductivity. 
Section VI is devoted to discussion. The mathematical 
details are presented in the appendices. 



A. Summary of the key results 

The most important results of Section |TT] are the ir- 
reducible representation of the space group PA/nmm at 
the F- and M-points. These are summarized in Tables U 
and mil The rest of the paper builds on these results. 

Our key results for the low-energy effective theory are 
the definition of the low energy effective "spinor" ip in Eq. 
(J17p whose components transform under the irreducible 



representations discussed in SectionlTTl Additionally, Eqs. 
(Ulll-dTni) (no spin-orbit), and Eqs. (EUl-dSTl) (with spin- 
orbit), correspond to the continuum description — con- 
sistent with the underlying crystalline symmetry — of 
the electronic states which are most affected by either 
magnetic or superconducting ordering. In the same sec- 
tion, the contact four-fcrmion interaction term is given 
in Eqs. pil)-(|5^. 

Section IIVI presents the symmetry analysis of spin 
density wave order parameters together with the elec- 
tronic spectrum and the Fermi surfaces for the collinear 
spin-density wave. Fig. I13[ and co-planar four-fold spin- 
density wave, Fig.Illl In Subsection IIV Al the conditions 
for the Kramers degeneracy are analyzed. In Subsec- 
tion |IVD] we provide the generic, symmetry based, argu- 
ments for the presence or absence of degeneracies in the 
spectrum using a generalization of Wigner-von Neumann 
analysis. 

The main finding in section |V] is that an s-wave super- 
conducting state can be well described with only three 
pairing terms, Eqs. ([5^ and (|57p . Although these terms 
arc k-independcnt, we find that while the pairing gap 
on the hole Fermi surfaces is isotropic, Eq. ([56|. it is 
k-dependent and anisotropic on the electron Fermi sur- 
faces with the possibility of gap nodes, Eq. ([M]). With 
the spin-orbit coupling, two additional, k-independent, 
spin triplet pairing terms, Eqs. ([73]) and (|80)) . are al- 
lowed. These lead to gap anisotropy on the hole Fermi 
surfaces, as seen in Eq. ([75]). The Eq. ([M|) gives the gap 
on electron Fermi surfaces in the presence of spin-orbit 
coupling and the concomitant spin triplet admixture. 



II. IRREDUCIBLE REPRESENTATION OF THE 
SPACE GROUP PA/nmm 

Our main analytical tool in this paper is the space 
group symmetry together with its irreducible representa- 
tions. We use the method of C. Herring^ to construct 
the irreducible representations of (non-symmorphio^) 
PA/nmm at the T- and M-points, and provide the key 
steps in what follows. The results thus obtained are first 
used to study the symmetry properties of the physical 
states and operators in iron-pnictidcs; they are then used 
to construct the low-energy effective model, the order pa- 
rameters, and to study the physical consequences of spin- 
orbit coupling on both the spin-density wave state and 
the superconducting state. 

Throughout this paper, we use the Seitz notation for 
symmetry operations, {3|t}, where g is a point group 
operation (rotation, reflection, etc.), which keeps the co- 
ordinate center invariant, followed by a translation by 
vector T. The product rule is 



integer) translations. 



{.91 1 Ti} {.92 1 7-2} = {5iff2|Ti +.91T2}, 



(1) 



T — \ {e|mim2TO3} mi, m2, 777,3 e Z j- 



(2) 



is an invariant subgroup of the space group, i.e., if we 
take any clement {glr} of PA/nmm, then for any 77ii's, 
the product {g\T}{e\mim2mz}{g\T}~^ is also an integer 
translation. The three numbers 777,^ represent a transla- 
tion by t = Y.i^i 'm-iO-i- 

PA/nmm is a non-symmorphic group. Therefore, it 
is impossible to choose a coordinate center such that 
every symmetry operation can be taken to be a point 
group operation followed by an integer translation, i.e., as 
{5|7Tii77i277i3} with all mi integers. This is because of the 
presence of the r7 glide plane which involves a fractional 
(non-integer) translations by tq = (22*^)' combined with 
the a6-plane mirror. 
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where 51X2 is the result of the point group element gi 
acting on vector T2. The group of all full unit cell (or 



FIG. 1: Bottom left: the coordinate systems used in this 
paper. Right: the lattice of iron based superconductors, rep- 
resented here by a single layer. Iron atoms (dark red) are split 
into two sublattices, A and B. Pnictide atoms (green) sit at 
iron-plaquette centers; one sublattice of the pnictide atoms is 
puckered above, and the other sublattice is puckered below 
the layer. One unit cell, centered on the coordinate origin 
(crosshairs), is outlined by the black square. The three gen- 
erators of Pr are shown in blue. The vertical mirror {(t^|00} 
passes through the coordinate origin. The vertical mirror a^ 
also passes through the coordinate origin, however, the sym- 
metry of the lattice requires this mirror to be followed by 
a translation by tq. The combined operation, {o"'^|^|}, is 
equivalent to a mirror which passes through pnictide atoms. 
The third generator is an n-glide mirror, {o"^|^^} and it has 
no fixed points. The action of each mirror is illustrated on 
an 'R' symbol. Under the operation of any of the first two 
generators, a mirror image of 'R' is created. The third gen- 
erator, acting on the 'R' above the plane, glides this symbol 
and puts it below the plane, which we represent by a hollow 
'R'. Top left: the first Brillouin zone for this lattice. 



There are 8 point group operations, {gjOOO}, centered 
on an iron atom (see Fig. [1]), forming the group Ti2d- 
The group can be generated by two elements: {cr^jOOO}, 



a mirror reflection about the yz-plane, and {5*41000}, a 
90° rotation about the z axis, followed by the mirror re- 
flection about the xy-planc. In addition to the eight point 
group operations, we also have an inversion followed by 
the fractional translation, {i|iiO}. Together they can be 
used to generate all the elements of PA/nmm. 

When we combine the 8 elements of D2d, with the 
8 symmetry operations obtained by multiplying {i|^^} 
with the elements of D2d, we find that this set is not 
closed under multiplication, and therefore, as is, it cannot 
form a group. (From now on, to simplify the notation, we 
omit the translation along the c-axis since it is through- 
out the paper). However, if we think of the product as 
being defined modulo integer translations, then these 16 
operations form a group, which we denote by Pp. More 
precisely^, Pp is the factor group {PA/nmm)/T. 

In order to find the irreducible representations of the 
space group, we first note that the irreducible repre- 
sentations for the Abelian subgroup of integer trans- 
lations, T, are specified by wave-vector k, such that 
-Dk({e|t}) = e'*'*. As usual, two irreducible represen- 
tations. Z^ki and I?k2j are equivalent if vectors ki and 
k2 differ by a reciprocal lattice vector, hence k's can be 
taken only in the first Brillouin zone. 

Since any element of PA/nmm can be casted as a prod- 
uct of an integer translation and an element in Pr, it is 
sufficient to determine how an object (e.g., a Bloch state) 
with momentum k transforms under all the elements of 
Pp. For each k, we first define the star of k as the set of 
all distinct momenta gk, where {^It} arc all the elements 
of Pr . The construction of irreducible representations of 
the space group is different when k lies inside or on the 
border of the Brillouin zone. 

When k is inside the Brillouin zone, the little co-group 
of wave- vector k, which we denote a Pk, is defined as the 
subgroup of Pr which keeps vector k invariant up to a 
reciprocal lattice vector. The irreducible representations 
of the space group arc then labeled by wave- vector k and 
another label for the irreducible representation of Pk. 
Any element of the space group which is a product of an 
integer translation and an element of Pk is represented 
as 

D^,,{{g\T = t + T'}) = e^'^*i?k. ({.gk'}) 

= e^'^-A({5|r'}), (3) 

where Di is a representation of Pk, and t is an integer 
translation vector chosen such that {<?|t'} is an element 
of Pk. Therefore, in order to determine symmetry prop- 
erties of a Bloch state with momentum k. it is sufficient 
to determine how it transforms under the members of the 
little co-group Pk, i.e., which irreducible representation 
of Pk it belongs to. Notice that the symmetry operations 
which do not keep k invariant change the momentum of 
a Bloch state from k to some k' which is in the star of k. 
Hence an irreducible representation at k also determines 
the symmetry properties for states with k' in the star 
of k, and once we find all the irreducible representations 
at some k, we have automatically found the irreducible 



representations at all other wave-vectors that belong to 
the star of k. (This is also true when k is at the border 
of the Brillouin zone.) 

If the momentum k labeling an irreducible representa- 
tion of the space group lies at the border of the Brillouin 
zone, and the space group is symmorphic, then the con- 
struction of the irreducible representations follows the 
same steps as when k is inside the Brillouin zone. 

When the momentum k is at the border of the Brillouin 
zone, and the space group is non-symmorphic, the irre- 
ducible representations of the space group at such a wave- 
vector cannot be constructed as described previously due 
to the presence of fractional translations. To demon- 
strate this, consider two symmetry operations, {gi|Ti} 
and {g2|'''2} and an irreducible representation Dk.i- The 
representation of the product of the two symmetry oper- 
ations, according to Eq. (jSj, is 

-Dk,i {{9i\ri}{g2\T2}) = -Dk.i ({5i52|n -f 31X2}) 

= e'^<^-+o^^-^D, {{gm\Ti + 91T2}) . (4) 

On the other hand, if wc take the representations for each 
symmetry operation and multiply them, we get 



i?M({5i|ri})i?k,a{52|r2}) 



zk-xi 



A({<?i|Ti})e^'^-^A({,92|r2}) 



e'''("^+"^)A({5i52ki+.giT2}). 



(5) 



These two expressions, Eqs. ^ and ([5|), must be equal, 
which is true if 

gik-(-ri +31-7-2) _ gik-(-ri+T-2) ^ gi-r2-(k-(gi)"ik) _ -^ iq\ 

for any gi and T2 . This is the case when the space group is 
symmorphic (r2 = 0), or k lies inside the Brillouin zone 
(then (gi)^^k = k). However, when the space group 
is non-symmorphic, and the wave-vector k sits on the 
Brillouin zone border, then {gi)^^\i and k can differ by a 
reciprocal lattice vector, in which case Eq. © is violated 
for some T2's. 

The irreducible representations of PA/nmm for k's on 
the Brillouin zone boundary are therefore constructed 
differently^. We discuss the irreducible representations 
of the space group at wave- vector M in Subsection IIIBl 
while the full construction of these is delegated to Ap- 
pendix |B] The irreducible representations at the other 
k-points sitting the Brillouin zone edge are enumerated 
in Subsection III CI 



A. Group Pr and its irreducible representations 

The 16 elements of the group Pr have been introduced 
earlier in this section. To repeat, 8 elements are given 
by {g|00} where g S T)2d, the other 8 elements are ob- 
tained from these by {i|^^}{(7|00} = {ig|^^}. These 
16 elements form a closed group Pp, where the Scitz 
multiplication rule, Eq. ([T|), is defined modulo an integer 
translation. 



The group Pr is isomorphic to D4/1. The mapping 
from Pr onto D4/J is performed by 'stripping-off' the 
translation part from the symmetry operation, {g|T} — )• 
g. Under the inverse mapping, an element g £ Ti^h 
is mapped onto {g|00} if g G D2d, and onto {gl^l} 
otherwise. We use this isomorphism when providing 
the irreducible representation tables, which arc widely 
available^, in this subsection. 



Pr 


{<^'\hh} 


i-'Wh) 


K|oo} 


f(k) 


^19/u 

^2g/u 
Blg/u 
B2g/u 

Eg/u 


±1 

Tl 
Tl 
±1 

[±1 1 
[ TlJ 


±1 
±1 
±1 
±1 

[ti 1 

[ TlJ 


±1 

Tl 
±1 

Tl 
r Tl' 

[ti _ 


g : k^ — ky 

u: {±k, + ky), 
i±kl + kl), 

rCxi^yyr^x ^ '^y ) 



TABLE I: The complete list of the irreducible representations 
at the F-point for the 3 point group generators. The addi- 
tional column on the right shows the symmetry properties of 
k-polynomials at kz — 0, up to the third order in k. Please 
note that the upper and the lower components of the axial 
vector representation Eg transform as Yz and —Xz, respec- 
tively. 

By the virtue of the isomorphism, we know that group 
Pr has three generators which, for future convenience, 
can be chosen to be {crxlH}, {cTzlii}, and {cr^jlOO}. 
The irreducible representations for the three generators 
are given in Table ID The irreducible representations for 
any other element of Pr can be obtained from these by 
multiplication. We use the names for the irreducible rep- 
resentations according to the isomorphism to D4ft,. The 
last column on the right in Table H] contains the classi- 
fications of the k-polynomials. Finally, Table |ll] shows 
the multiplication identities for the irreducible represen- 
tations of Pr- 



B. Group Pm and its irreducible representations 

The corner of the Brillouin zone, M, has a particular 
importance in iron based superconductors; the electron 
Fermi surfaces are centered around this point, as found 
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Aib 


Bit 


B2b 


Eb 


Ala 


Ale 


A2c 


Bic 


B2c 


E, 


Ala 
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Ai, 


B2c 


Bic 


E, 


Bla 


Blc 


B2c 


Ale 


A2c 


E, 


B2a 


B2c 


Bic 


A2c 


Ala 


E, 


Ea 


E, 


E, 


Ec 


Eo 


Ale © A2c © Blc © B2c 



TABLE II: The product table for the irreducible representa- 
tions of Pr — D4(i. The parity of the product is even, c = g, 
if the parities of the two multiplying irreducible representa- 
tions are the same, a = &; if the parities are opposite, a ^ b, 
the product has odd parity, c = u. 



thcorcticallj*iii2£i^"— and in experiments^"—. Follow- 
ing Eq. (J6|), we have demonstrated that the irreducible 
representations at the k-points at the Brillouin zone 
boundary, such as M, are not given by the little co-group 
Pk, a subgroup of Pp. Instead, there is a separate pro- 
cedure, due to C. Herring^; we outline the procedure 
here and give the final result: the list of the irreducible 
representations of PA/nmm at the wave-vector M. The 
technical parts of the construction arc shown in Appendix 

El 

The basic idea of the approach is to notice that the 
naively constructed little co-group at M, which has 16 
elements of Pr is not closed under the multiplication as 
demonstrated in Eq. (jS]). Therefore, the group Pm must 
contain additional elements. In this case, wc construct 
group Pm as a quotient group of PA/nmm with the in- 
variant subgroup of even translations 



Tm = {{e|t} cxp(iM-t) = l|. 



(7) 



Group Pm defined in this manner is closed under the 
Seitz product Eq. ^ modulo even translations. It con- 
tains 16 elements originally found in Pr, and additional 
16 elements obtained by a multiplication with an odd 
translation, {e|10}{5|T} == {^llO + r}. 

Any element of PA/nmm can be written as a product 
of an integer translation and an element of Pm • Since all 
even translations arc represented by the unity per defi- 
nition Eq. (I7|), it is sufficient to find the irreducible rep- 
resentations of Pm in order to determine the representa- 
tion for any element of PA/nmm. There is, however, one 
subtlety: the odd translation, {e|10}, which is among the 
elements of Pm , rnust be represented by 



DM4{e\lO}) 



iaWl-St-r 



-1. 



(8) 



Therefore, only the irreducible representations of Pm for 
which Eq. ([8]) is satisfied are physical at the M-point. 
All other irreducible representations are unphysical at 
M and should be disregarded. 

The group Pm constructed here is not isomorphic 
to any three-dimensional point groups, all of which 
are listed in Ref. |4^. We are therefore left with the 
task of deriving the irreducible representations of Pm- 
The full procedure, based on the method of induced 
representations^, is spelled out in Appendix |B] Here we 
list only the physical irreducible representations of Pm 
in Table IIIII The complete list of the irreducible repre- 
sentations is found in Appendix IB] 

All the physical irreducible representations at the M- 
point are two dimensional; the two components of each 
doublet are related to each other by the mirror reflec- 
tion {cr^ 1 00}. One immediate consequence is that all 
the states with momentum M must be doubly degen- 
erate; other consequences of this double degeneracy are 
discussed in the remainder of this subsection. 

The product table for the physical irreducible repre- 
sentations of Pm is given in Table IIVI A product of 



-Emi 

-Bm2 
E-M.3 
E-M.4 



Wll^} 



-1 

-1 
■ 1 0' 

1 

1 o" 

-1 

-1 

1 



{--1^^} 



-1 

1 
-1 

1 
-1 

1 
-1 

1 



K|oo} 



1 

1 

1 

1 

1 

1 

1 

1 



rj-taTTrpD 



COS / sin 
sin / cos 
sin / sin 
cos / cos 



TABLE III: The list of physical irreducible representations at 
the M-point for three group generators. Since PA/nmm is 
non-symmorphic, all the physical irreducible representations 
are two dimensional. This implies the double degeneracy of 
any states at the M-point which are classified according to 
these irreducible representations. The last column in the ta- 
ble should be used in Eq. (|10|l in order to construct the com- 
plete symmetry adapted basis of r-functions which transform 
according to one of the E^a representations. The illustrative 
examples of such functions are given in Figs. [5]and|31 
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TABLE IV: Product table for the physical irreducible repre- 
sentations of Pm. Since two vectors M add to zero in the 
reciprocal space, the product has no momentum, hence it de- 
composes onto irreducible representations of Pr — 04;,. 



two Bloch states, each with momentum M, has no to- 
tal lattice momentum, hence the products of two -Em 
irreducible representations are decomposed into the irre- 
ducible representations of Pp. The mixed product table, 
where one irreducible representation is of Pr , while the 
other is a physical irreducible representation of Pm, has 
the total momentum M. Therefore, such a product de- 
composes into -Em irreducible representations, as shown 
in Table El 

We use the physical irreducible representations of Pm 
for the symmetry classification of two-dimensional Bloch 
functions with momentum M, i.e., all the functions /(r), 
where r = (x, y, 0), such that 



/(r) = e^^-u(r) 

_ jMt V u e 



i=^{inix-\-rn2y) 



mi,m2G'Z 



E 



oi^(nia;+Ti2i 



(9) 



ni,n2eZ+l/2 



Here, u(r) is a lattice periodic function, with Fourier har- 
monics coefficients Wmi,m2- The function is alternatively 
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TABLE V: Product table between the irreducible representa- 
tions of D4;j and the physical irreducible representations of 
Pm. 



decomposed into the half-integer Fourier harmonics in 
the last line of Eq. (jH]). The half- integer Fourier harmon- 
ics form a complete basis for the two-dimensional M, 
however, they are not adapted to the lattice symmetries. 
To remedy this, we use the projector method on the half- 
integer Fourier harmonics, and find a new basis such that 
each basis function transforms as a particular component 
of one of the -Em's. All the functions in this symmetry 
adapted basis come in doublets, as they should, due to 
the two-dimensionality of each Em • For any physical ir- 
reducible representation, the two doublet states can be 
written as 
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(10) 



where integers 7711.2 > 0, and T^ '" are two trigonometric 
functions, listed in Table ITTll for each physical irreducible 
representation. The basis defined by Eq. ([TU]) is complete. 

The lowest harmonics, mi = 1712 = 0, are non-zero 
only for Em2 an Em4- These four functions are plotted 
in Fig. [21 In the same plots, iron atoms are marked by 
solid dots, and one unit cell is outlined as a guide to the 
eye. For any higher harmonics (i.e., 7711+7712 > 0), and for 
each Em, there is precisely one doublet of functions. For 
illustration, we plot the second lowest harmonics (r^i = 
0, 7772 = 1) in Fig. O 

In both Figs. |2l and |3l the Em4 harmonics have minima 
and maxima on iron atoms and nodes on (the a-b plane 
projection of) the pnictidc atoms. All other harmonics 
vanish precisely on the iron atom sites. Similarly, the 
Em2 harmonics have the largest magnitude on pnictide 
atom sites, while all the other harmonics vanish there. 
This is true for any higher harmonics. It can be demon- 
strated from Eq. ^TU\i that only Em4 harmonics are node- 
less on iron-sites; on pnictide sites only Em2 harmonics 
have no nodes. Experimentally, the spin-wave density 
order is found to have magnetic moments located on iron 
We conclude therefore that the orbital part 



atom; 



,32.33 




FIG. 2: The lowest order harmonics (mi = m2 = 0) in the 
symmetry adapted basis of functions given by Eq. (|10p . No- 
tice how the E-M4 functions have a finite value at the positions 
of iron atoms and vanish on pnictide sites. Conversely, the 
Em.2 functions are finite on pnictide atoms and vanish on iron 
atoms. 





M2 J^M3 




ing layers^i^. Therefore, we should consider the sym- 
metry adapted fmictions at A = 7r(l/a, 1/a, 1/c) instead 
when discussing the spin-density wave order found in the 
experiments. Since A is at the edge of the Brillouin zone, 
Pa must be constructed according to the same procedure 
as Pm- It turns out that these two groups are isomor- 
phic. Pa — Pm; and, as a consequence, their physical 
irreducible representations must be equivalent. There- 
fore, the symmetry properties of spin-density waves in 
iron based superconductors are, strictly speaking, gov- 
erned by i?A4, which is given by the same matrices as 
Em4 in Table UTTl 



C. Irreducible representations of P4:/nmm at 
lower-symnietry points in the Brillouin zone 

In this subsection we focus on the remaining symmetry 
indistinguishable k points in the (two-dimensional) Bril- 
louin zone. These are X = {tt/q, 0, 0), lines S = F — M, 
A = F - X, and Y = X - M, and the set of all other 
k-points not sitting on any of the high-symmetry points 
or lines. The expansions around the F- and M-points, 
discussed in the previous two subsections, are, in prin- 
ciple, independent of the classification provided in this 
section. Nevertheless, the results presented here will be 
used in comparing our results to two- and three-orbital 
tight-binding models, as well as, in understanding the 
appearance of nodal points in the symmetric, and in the 
coUinear spin-density wave state. 

The high symmetry point X = (tt/a, 0, 0) lies at the 
edge of the Brillouin zone. Since PA/nmm is a non- 
symmorphic group, the point group at X is constructed 
in a similar fashion as that at M. We do not give the 
details of the construction in the paper; instead, we only 
present the final result: group Px is the subgroup of 
Pm whose elements keep X invariant. This group has 
16 elements and is isomorphic to Jy^h- An isomorphism 
between the two groups is defined through the generators. 



{fj^lOO} 

{^2122} 
{a^|10} 



X 



(11) 



FIG. 3: For each physical irreducible representation E-mz, we 
plot the symmetry adapted functions doublet. In this figure 
we show the next lowest order harmonics (mi — 0, m2 = 1) 
given in Eq. HlOf) . Similarly to Fig. [21 the only doublet with 
non- vanishing values on the iron sites is -Bm4; the only doublet 
with non- vanishing values on the pnictide sites is E-M2- 



of any spin- wave density wave order in iron based super- 
conductors must have the symmetry properties of i?M4- 

The last statement comes with one caveat. In the 
present exposition, we constrain ourselves to the analy- 
sis of a single layer. The experiments also show that the 
spin-density wave has the opposite sign in two neighbor- 



All the other elements of Px can be generated from these; 
as a consequence of this isomorphism {e|10} — >■ C2- Just 
as it was the case with the Pm group, the only physi- 
cal irreducible representations of Px are those for which 
D({e|10}) = -1. Using Table U we find that there 
are only two such representations, -Exg and -Exm- Since 
both physical irreducible representations of Px are two- 
dimensional, one consequence is that all the bands at this 
point must be doubly degenerate. To distinguish between 
the two irreducible representations, one easy method is 
to observe the parity of the Bloch states at X under the 
mirror reflection {cr^|00}. The Bloch states transforming 
according to i?xg are even, those transforming according 
to E'xm arc odd under this operation. 



The group Py is the subgroup of Px which keeps a 
wave- vector k G y invariant. This group is isomorphic 
to Civ, with the isomorphism defined by the first two 
Hues in Eq. (|lip . The mirror reflection {ct^|00} is ab- 
sent in Py. The F-hne constitutes the Brillouin zone 
boundary, and there is only one physical irreducible rep- 
resentation, E, specified in Table IVll The symmetry clas- 
sification at the K-line is therefore almost trivial since 
everything transforms according to the same irreducible 
representations. However, one should notice that, since 
E is two-dimensional, each doublet of Bloch functions at 
\s. G Y must be chosen in such a way that the two dou- 
blet components transform into each other (up to a phase 
factor) under the {<t^|||} mirror reflection. One conse- 
quence of this double degeneracy is that all the Bloch 
states at any k £ Y arc doubly degenerate. 



Py 


a"" 
{a" 100} 
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-1 
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-1 
-1 



TABLE VI: The unique physical irreducible representation of 
the generators of point group Py = C4„ . 



The remaining k-points are located inside the Brillouin 
zone. For any of these, group Pk is the little co-group of 
Pr- The representations of an element of PA/mnm are 
given by Eq. ([3]). 

For k € S, the little co-group Ps is generated by two 
mirror reflections, {cr^lH} and {(T^|||}. This group is 
isomorphic to C2i,, with the isomorphism 



KIHI 



KIH} 
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The table of the irreducible representations and their 
product table are given in Table FVlIl 

At last, for a state with a momentum k = (fc^;, fcj,, 0), 
lying on none of these high symmetry points or lines, 
there is still a single mirror reflection {c^looi which 
keeps k invariant. This is due to the fact that we con- 
tained ourselves to the two-dimensional Brillouin zone 
and kz = 0. This mirror defines group Pk = Cg for a 
general k-point. Group Cg has two irreducible represen- 
tations A' and A" corresponding to the states that are 
respectively even and odd under the action of {c^lHj 
mirror reflection, in sense of Eq. ([3]) . 

This list exhausts all the possible k-points in the two- 
dimensional Brillouin zone for iron-pnictides where kz = 
0. ft appears that we omitted the irreducible representa- 
tions of the PA/nmm with wave- vector k = {kx, ky, tt/c). 
These irreducible representations are relevant for the 
Bloch states or order parameters which change the sign 
in alternating layers, as mentioned in Subsection III Bl It 
turns out that, for the PA/nmm group, the irreducible 
representations at any k = (kx, ky, tt/c) are equivalent to 
those at k = [kx, ky,0). Therefore, the irreducible repre- 
sentations at the F-point are also governing the symme- 
try properties at Z = (0,0, tt/c), the physical irreducible 
representations at R = (tt/a, 0, tt/c) are the the same as 
those at the X-point, and most importantly, the physi- 
cal irreducible representations at A = (tt/o, 7r/a, tt/c) are 
equivalent to the four E^s. Similarly, the three high- 
symmetry lines S, A, and Y, share the symmetry clas- 
sification with their kz = tt/c counterparts, U-, S-, and 
T-line, respectively. For a general k ~ {kx,ky,TT/c) mo- 
mentum, Pk = Cs- Although kz — tt/c implies that k 
lies on the Brillouin zone border, there is no issue with 
the non-symmorphicity of PA/nmm here since the frac- 
tional translation vector Tq has no z-component. 



The table of the irreducible representations and their 
product table are given in Table IVIII 
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TABLE VII: The irreducible representations for the gener- 
ators of the isomorphic point groups C2u = Ps — Pa- 
These are used in the classification of Bloch states with k 
atE = r-MorA = r-X lines in the Brillouin zone. The 
product table is given on the right. 

For k € A, the little co-group Pa is generated by two 
mirror reflections, {tr'^lOO} and {ct^|||}. This group is 
isomorphic to C2i,, with the isomorphism 



K|00} ^ fj^ 



K|ii}^a^ 
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D. The classification of the Bloch states near Fermi 

level 



The first step in constructing an effective theory is rec- 
ognizing which states, in the full model, are the most 
physically relevant and which states bring only quanti- 
tative corrections and may be integrated out. Gener- 
ally, states with the energy at, and in the vicinity of, the 
Fermi level are the ones to be kept in the effective model. 
The dispersion of these states should match the disper- 
sion of the states of the full model as closely as possible, 
but without overburdening the effective model with un- 
necessary details. Perhaps even more importantly, the 
states of the effective model must have exactly the same 
symmetry properties as the states in the full model they 
originate from. Otherwise, the model may not be able 
to capture the correct symmetries of physical effects in- 
tended to be described by the effective model. For ex- 
ample, the expectation value of a bilinear may not be 
associated with an order parameter of correct symmetry, 
transitions violating selection rules may occur, etc. If we 



are going to construct an effective theory for iron based 
superconductors, we ought to identify the symmetries of 
every state that will make a part of the theory. That is 
the main result of this subsection. 

The symmetries of states can, in principle, be deter- 
mined from a 'first principle' band structure calculation. 
A direct output of such a calculation is the composition of 
eigenstates and from there, the symmetry nature of each 
eigenstate may be derived. An alternative, and the one 
that we use here, is the band structure obtained from a 
tight-binding model. We use two independently derived 
tight-bindin g rn odels. The details of the first model are 
given in Ref. ^^ That tight-binding model uses 16 states 
per unit cell, five d states per each iron atom and three p 
states per each pnictide atom, and is fitted to the band 
structure of undoped LaOFcP obtained in an ab initio 
calculation^!. The other tight-binding model we use is 
presented in Ref. |42|. That tight-binding model is fitted 
to a first principle band structure calculation for a; = 0.1 
doped LaOFeAs. The band structure calculation is pre- 
sented in the same paper. This tight-binding model uses 
five maximally localized Wannier states per iron atom, 
each of these states having orbital symmetries of one of 
the iron 3d orbitals. Between the two tight-binding mod- 
els, we find no difference in the symmetries of the bands 
which cross the Fermi level. 

We use the irreducible representation of PA/nmm, pre- 
sented in this section, to classify the iron states at any 
k. In a tight-binding model, an iron state at momentum 
k is given by 



always, split into two one dimensional irreducible repre- 
sentations once the symmetry is lowered. 



^k(r)=EE' 



,Jk-(R,-(-5) 



^0^(r-R,-(5)4(k).(14) 



Here, R; are the positions of iron atoms in the A sub- 
lattice; 5 takes two values, or tq corresponding to the 
A or B lattice, respectively, ^^^^(r) is the orbital part of 
the wave function for fj., an iron 3d orbital centered at 
r = 0. There are ten iron states for each k, two per iron 
atom in a unit cell, and their symmetry properties are 
presented in Table rvTIIl 

These tables contains additional information related to 
the symmetry lowering in the vicinity of high symmetry 
point or lines. Namely, if a band is classified according to 
a certain irreducible representation at a high symmetry 
point, then its symmetry properties in the vicinity of that 
point can be deduced from there. This is a consequence of 
the fact that the transformation properties of that band 
under the remaining symmetry operations, away from the 
high symmetry point, arc inherited from the properties 
at the high symmetry point. Conversely, the symmetry 
properties at the high symmetry point may be seen as 
the same properties from its vicinity augmented by addi- 
tional symmetries at the high symmetry point. The same 
is true when the symmetry is lowered from a high sym- 
metry line to its vicinity. One dimensional irreducible 
representations at a high symmetry point uniquely de- 
termine the symmetry properties in its vicinity. The two 
dimensional irreducible representations usually, but not 
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TABLE VIII: The symmetry classification of the tight-binding 
states, Eq. (|14[) at the three high symmetry points, F, M, and 
X, is shown in the second column of each table. In the vicinity 
of each high symmetry point, the symmetry is lower than it 
is at the high symmetry point. The remaining columns show 
the symmetry properties of the same states inherited from the 
symmetry properties (i.e., the irreducible representation) at 
the high symmetry point. Two-dimensional irreducible rep- 
resentations usually get split into two one-dimensional ones. 
Physically, that corresponds to an opening of a gap due to a 
reduced symmetry. 

Having recognized the symmetry properties of each 
tight-binding iron state, we diagonalize the Hamiltonians 
in Refs. [ll| andl43 and establish the symmetry properties 
of each eigenstate based on its mixture of iron orbitals. 
We have performed this classification for any k in the 
Brillouin zone. For clarity, in Fig. 2] we present the band 
structure along the F — M (the S-line) . As shown in the 
legend of the plot, the color/dashing coding for each band 
corresponds to the irreducible representation of Ps which 
gives the symmetry properties of that band. On the left- 
and right-hand side of the plot, the irreducible represen- 
tations for the bands at F- and M-points, respectively, 
are shown. Notice that the irreducible representations 
for each band at the high symmetry points and along the 
S-line are matching in accordance with Table IVHII The 
Fermi surfaces for these two tight-binding models and 
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the symmetry properties of the Fermi surface states are 
shown in Fig. [5] The bands and their symmetry prop- 
erties shown in Figs. H^ and Hb have the same outhne. 
The major difference is the order of bands at the F-point; 
however, it does not concern the bands which cross the 
Fermi leveL These bands have the same properties not 
only along the E-line, but at any momentum k, as seen 
from Fig. \5\ 
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form a degenerate E'g-doublet. Along the high symme- 
try lines, S and A, the doublet splits into one A2 and 
one B2 state. These two states are not degenerate un- 
less the Hamiltonian is fine-tuned. The upper band, re- 
sponsible for the outer Fermi surface, has the symmetry 
properties governed by A2 irreducible representation. It 
is odd under {c^l^^} mirror, and also odd under {cr^|00} 
or {cr'^l^^} mirror. The other band, defining the inner 
hole Fermi surface, has symmetry properties given by 
B2 irreducible representation, hence it is odd under the 
{(T^|i^} mirror, and even under the other mirror of Ps 
or Pa group. Away from the high symmetry lines, the 
states on the two hole bands are classified as A" , i.e., odd 
under {ct^I^^} mirror, neglecting the phase due to the 
translation part, Eq. ^. 





A" 


m 


A2i^) 


m 


S2(A) 


• 


A2i^) 


• 


^2(2) 





-- A' 
A" 

• ^2(2) 

■ E(Y) 



r,ig- 
Aig- 
Bzu' 



. Jv 



'^ 



°2gt 



Byil 

~' -Bm3 



~" ■■£m4 
-Bm2 



M 



b) 



FIG. 4: The band structure along the E = F — M line in 
the Brillouin zone (see Fi g. P reconstructed from the tight- 
binding model in Refs. a) [ll|, b) 142 . On the left-hand side of 
the plot, the symmetry properties of the bands at the F-point 
are shown; analogously, on the right-hand side, the symme- 
try properties for each band at the M-point are shown. In 
between, the color/dashing corresponds to the symmetry clas- 
sification of the bands according to Pe — C21, group. 

The symmetry properties of each band crossing the 
Fermi level in iron-pnictides are now well understood. 
The two hole bands originate at the F-point where they 



FIG. 5: Fermi surfaces, reconstructed from the tight-binding 
model in Ref. |Tl|. The color of each Fermi surface at an arbi- 
trary momentum k corresponds to the parity of states at the 
surface under the {o"^|||} mirror. Along the high symmetry 
directions the symmetry is augmented, and the states on the 
Fermi surfaces are classified accordingly. 



The symmetry properties of the electron bands in the 
vicinity of the M-point are quite different. Let us first no- 
tice that the two electron bands in the vicinity of the M- 
point originate from the £'mi and E'ms doublet. Along 
the S-line in the Brillouin zone the doublet with the 
higher energy at M, £'mi, splits into one A2 and one B2 
state. The A2 state is weakly dispersive and, although 
it connects to the Eg doublet at the F-point where it 
crosses the Fermi level, we can consider that this state 
does not cross the Fermi level in the vicinity of the M- 
point. The other state, i?i, has a steep dispersion such 
that it crosses the Fermi level close to the M-point. This 
is the state which corresponds to the point on the minor 
axis of one of the electron Fermi surfaces. Observing the 
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other doublet, £'m3, we find that along the E-line it splits 
into two states, Bi and i?2- The Bi state originating in 
the -Ems doublet has a nearly flat dispersion and it does 
not cross the Fermi level. The other state, B2, does cross 
the Fermi level, and this state corresponds to the point 
on the major axis on the second elliptical Fermi surface 
near the M-point. As one moves away from the high sym- 
metry line S, the symmetry is reduced to C^. The two 
Bi states arc both even under a^ mirror, i.e., belong to 
the A' irreducible representation. The two other states, 
A2 and B2, are odd, i.e., transform according to A2 ir- 
reducible representation. Notice that these two states 
therefore mix at any k which does not lie on the S-line. 
However, precisely at the E-line, these two states split 
into two different irreducible representation and their re- 
spective bands must cross by the symmetry constraints. 
Hence, these two bands form a Dirac point on the S-linc. 



By extending the analysis of the states to the edge of 
the Brillouin zone, the y-line, we notice that the two up- 
per bands form one i?-doublet that originates from the 
i^Mi-doublet at the M-point. It is the states of this dou- 
blet that form two degenerate points where the Fermi 
surfaces intersect. One state is even and the other odd 
under the cr^-mirror. Similarly, the two lower bands, the 
ones that originate from the i?M3-doublet and do not 
cross the Fermi level, form an £'-doublet along the Y- 
line. Crossing into the second Brillouin zone, we wish 
to connect the degenerate states from the F-line to the 
states on M — F'-line, where F' = (27r/a, 0). It is im- 
portant to notice that the splitting of the i?M-doublets 
changes along this line as compared to the E-line. The 
reason is that along the M — F'-line, the two symme- 
try mirrors are {o"^|ii} and {c^j^^}. Compare that to 
the E-line where the first mirror was {o'^|2 2}- Conse- 
quently, the E'Mi-doublet splits into one Ai and one B2 
state. The B2 state is strongly dispersive and it crosses 
the Fermi level close to the M-point. The Ai state is 
weakly dispersive and it does not cross the Fermi level in 
the vicinity of the M-point. Even though the symmetry 
generating mirror is different along the M — F'-line, the 
i?M doublet is still split into one Bi and one B2 state. 
The difference, with respect to the E-line, is that the Bi 
state emanating from the -Ems crosses the Fermi level 
while the B2 state is weakly dispersive and it does not 
cross the Fermi level. 



III. LOW ENERGY EFFECTIVE THEORY AND 
SPIN-ORBIT COUPLING 

In the previous section we have identified the sym- 
metry (i.e., the irreducible representation) of the states 
within ~ 2eV of the Fermi level. This information, and 
the analogous classification of k-polynomials, allows us to 
form all invariants that can appear in an effective Hamil- 
tonian. We now construct such a low energy effective 
Hamiltonian accurately describing the states which cross 
the Fermi level. 

The non- interacting part of the Hamiltonian in the nor- 
mal state is 



'Ho — Hq 



H. 



(15) 



The second term accounts for the spin-orbit coupling, 
and we discuss it in the next subsection. The first term 
is spin 5(7(2) symmetric, and reads 



Ho 






hU^) 




M^), 



where the six component (pseudo) spinor is 



V'.(k) 



^x,.(k) \ 
V'y,a(k) . 

#,a(k) ; 



(16) 



(17) 



For each spin projection, the upper component of ipx 
transforms as E^ and the lower as E^ ; the upper 
component of tpy transforms as E^ and the lower as 
E^ . Similarly, for each spin projection, the upper and 
the lower components of the spinor '0r(k) transform un- 
der the (axial vector) Eg representation at the F point as 
Yz and —Xz, respectively. The 2x2 blocks in Eq. ITF 
are 



hU^) 



/ir(k) = 



iv±(k.) 



f-r + ^ h bkxky 



-iw±(k) 



es + ^iag/c^fcy 



c kl 



■■{kl-kl) er + ^-bk,k 



(18) 



(19) 



The bands with states that belong to the Ai or Bi irre- 
ducible representations at the high symmetry line M— F', 
have all their states even under the {a^ I 2 5} niirror, hence 
away from this line the states transform according to the 
A' irreducible representation of Cs group. Conversely, 
the bands whose states at the M — F'-line transform ac- 
cording to the B2 irreducible representation are odd un- 
der the a^ mirror, and the states on these bands trans- 
form as A" at an arbitrary momentum k away from the 
high symmetry lines. 



with 



i'±(k) = V {±kx + ky 



■Pi 



±kt 



' P2^xf^y \^x ^ f^y) • 



(20) 



In the above, k is in-plane and measured in units of the 
inverse lattice spacing. The form of this Hamiltonian fol- 
lows from the symmetry properties of "0 and the product 
tables |TT] and IIVI for the irreducible representations at F 
and M. Finding the spectrum of Ho is now reduced to 
solving a simple 2x2 eigenproblem. 
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Note that we use the exact Bloch cigenstates at only 
two k-points, F and M. By the very nature of this k • p 
construction, we always deal with an analytic expansion 
of the effective Hamiltonian in powers of the components 
of k. This is not the case within other approaches to this 
problem^i^i^"— . There one first constructs the exact 
Bloch cigenstates at each k, and then uses such states as 
the basis for the expansion of the electron creation and 
annihilation operators. Among the many advantages of 
our approach is its suitability in a study of the effects 
of a quantizing magnetic field, since in Eqs. (|18p and 
([T^ the electro-magnetic vector potential A can now be 
minimally coupled^^. While straightforward in principle, 
exploration of such effects, however, is beyond the scope 
of this paper. 

To find the undetermined coefficients, wc first diago- 
nalize the tight-binding Hamiltonian at F and then use 
the resulting eigenvectors to construct the matrix ele- 
ments of the Hamiltonian at an arbitrary k away from F. 
The off-diagonal matrix elements vanish as k vanishes. 
Therefore, we can integrate out all the states except for 
the Eg doublet (see Fig. 2]) and expand^ the resulting 
2x2 effective Hamiltonian to second order in k. An anal- 
ogous procedure was carried out at M, except we kept 
the Emi doublet and the Em.3 doublet, making the effec- 
tive Hamiltonian at the M-point 4x4. Still, our choice 
of basis in Eq. (J17p , and our ignoring of any out-of-plane 
/c^-momentum dispersion, allows us to write such 4x4 
Hamiltonian in a block diagonal form. Depending on 
the tight-binding Hamiltonian we use, Ref. [Ill or Ref. |42, 
such method results in the values tabulated in Table HXl 
and stated in meV's. 
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TABLE IX: The parameters of the low energy effective model, 
Eq. p8ll9[) . derived from the two tight-binding models refer- 
enced in the text. 



b) 



FIG. 6: The comparison between the band structures for the 
tight-binding model (black, dashed) and the low energy effec- 
tive theory (red, solid) along the E-line in the Brillouin zone. 
The two panes correspond to a) Ref. [ll| and b) Ref. 14^. Com- 
paring the band structures in the vicinity of the Fermi level, 
we find that the low energy effective theory matches the tight- 
binding band structure almost perfectly near the F-point and 
the two are close in the vicinity of the M-point. 



The spectrum obtained from Eq. P^ (red, solid) 
is compared to the full tight-binding dispersion (black, 
dashed) along the E-line in Fig. [51 The Fermi surfaces 
are compared in Fig. [T] For the hole bands, the Fermi 
surfaces are captured with a marked precision; the hole 
bands dispersion of the low-energy effective model does 
not deviate significantly from the dispersions of the tight- 
binding models for the energies of order of few hundred 
meV's from the Fermi level. For the electron bands, the 
shape of their Fermi surfaces is well captured by the low- 
energy effective model, although there are small quanti- 
tative deviations. The dispersions too have small devia- 
tions from the tight-binding bands in the vicinity of the 



Fermi level. The fit may be improved by including higher 
order terms (quartic) in k in Ecj^s. (|T8|) and (fT9|) . 



A. Spin-orbit coupling 



In order to include the effects of spin-orbit coupling, we 
begin by analyzing the transformations of charge neutral 
Fermion bilinear operators under the symmetry opera- 
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FIG. 7: The comparison between the Fermi surfaces in the 
tight-binding model (black, dashed) and the low energy effec- 
tive theory (red, solid). The two panes correspond to a) Ref. 
[nl and b) Ref. \M- 



tions R = {.g|T}. We have 

(21) 



In order to find the terms in Hgo, i.e., the invariants 
of the form M^ ^ij;'^ a'^aij^^ where the repeated indices 
arc summed over, we look for the identity representation 
in the product of the three representations on the right 
hand side of Eq. ((2T|) . 

The axial vector representation of D4/1 is decomposed 
into A2gOEg. The z-componcnt of the axial vector trans- 
forms according to the A2g representation, while a^ and 
a^ transform as an Eg doublet. Recall that the low- 
energy states at the F-point transform as an Eg doublet 
of Dih- We can therefore form four independent bilin- 
ears composed out of these two states which transform 
according to Eg ^ Eg = Aig © ^23 © ^ig © ^2g- It fol- 
lows that there is only one k-indepcndcnt invariant in 
D"^^ ® Eg® Eg and that it only contains a^ . 

Similarly, for the low energy states at the M-point, we 
have i?Mi (8) E'mi = ^ig © ^23 © ^2« © -Bim and E-Ms <X) 
Em3 = Aig ® B2g © Alu © B2u- Therefore, there are 
no symmetry allowed, k- independent, spin-orbit terms 
with both Tp* and ijj transforming under £'mi , or both 
transforming under E-m^. On the other hand, ii^Mi ® 
Em.3 = Eg ® Eu- Therefore, there is a k- independent 
spin-orbit term mixing i^Mi and i^Ma states, containing 
the Pauli spin matrices a^ and tr^, but not cr^. 

Putting it all together, we find that 







uso 
''■M.aa' 



In the above, the sum over a =t, i and l3 =t,l should Hso — ^^ipl(k) [hf^^ j ^^^ 



be understood. Here, cr^^ is either a Pauli spin ma- 
trix ^ = j = x,y,z or a unity matrix, a^ ~ I2. 
In Eq. ()21|) a, b, c and d may be any one of the six 
components the spinor. For each R, the non-zero cl- 

D.iiR) Di3(i?) 



k cf..(t' 






uso 

"r,(T(T' 



V'<.'(k), 



(25) 



ements of matrix Dab{R) are 

E-M.i{R), 



D55{R) 

De^iR) 



D22{R) 

Di2{R) 
Dz^iR) 
DeeiR) 



D2i{R) 

D44{R) 



D^iiR) DssiR) 

= EM.3iR), 



and 



Eg (i?) . The transformation 



where 



and 



.. = ^Ar 



-i 
i 



properties of the Pauli spin matrices a^ follow from 
cr° = 1 corresponding to spin singlet, and a^, a^ and 
a^ corresponding to spin triplet. Naturally, a spin sin- 
glet is left invariant under all space group transforma- 
tions, hence for each R we have Dqq^ (R) = 1 and 
D^/{R) = Df^[R) = 0. On the other hand, the spin 
triplet operators transform according to the axial vec- 
tor representation of the point group component R in 
R — {.g|T}, and are left unaffected by the translation by 
T. Therefore, Dfy [R) — gf^ ■ The action of various 
symmetry operations in PA/nmm on the spin operators 
a^ can now be readily obtained from the action of the 
group generators. For /i = i and ly = j, 



'SO 



\ 
2AM 



1 






1 



(26) 



(27) 






D^^{{a^\00}) 



diag(l,-l,-l), 
diag(-l,-l,l). 




(22) 
(23) 

(24) 



Eq. ((26)) resembles the term introduced by Kane and Mele 
in the single layer graphene^. In fact, at the F-point 
there is a close similarity with the effective Hamiltonian 
for the bilayer graphene, the difference being the large 
particle-hole asymmetry term proportional to k^ which 
is negligible in the bilayer graphene. 

In addition to these, k-independent, spin-orbit terms, 
the symmetry of the lattice allows several other k- 
depcndent terms in the Hamiltonian. We may, however, 
neglect these because they originate from higher order 
hopping processes and are much smaller than the on-site 
terms, Eqs. dMI and (P7)) . 

If Figs. [5^1 andlHjD, we compare the dispersion obtained 
using the low energy effective Hamiltonian JIq with and 
without the spin-orbit term Hso in Eq. ([^5)1 . We use the 
values Ar = Am = 80meV obtained from from Ref. [2^. 
Both with and without spin-orbit coupling interaction. 
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we find that each band is doubly degenerate. This de- 
generacy is an obvious consequence of the spin double 
degeneracy in the absence of spin orbit coupling. When 
the spin-orbit coupling is present, the double degener- 
ate states form Kramers doublets — two states related 
to each other by a combined operation, time reversal Q 
followed by "glide" -inversion {j|^^}, must be degenerate 
since this operation is a symmetry of the Hamiltonian at 
any k-point in the Brillouin zone. Therefore, each band 
in Figs.|S^ and[5}D, as well as generally in iron-pnictides, 
is doubly degenerate. We find that, in the presence of 
spin-orbit interaction, any two bands never cross except 
at the M-point, where we find only Dirac-like points (see 
Fig. [5Jd). The Dirac points are guaranteed by the space 
group symmetry and cannot be removed without lower- 
ing this symmetry. 
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FIG. 8: The dispersion of the low-energy effective theory 
without (thin-dashed lines) and with (thick-solid line) spin- 
orbit effects included. The two panels correspond to the ef- 
fective theory near the a) F-, and b) M-point. The strength 
of the spin-orbit couphng we use is Ar = Am = 80meV. Each 
band is doubly degenerate regardless the presence of the spin- 
orbit coupling. No two bands cross except at the M-point 
where bands meet in pairs and each state is four-fold degen- 
erate. 



B. Comparison to other models for iron-pnictides 

After wc have constructed our low-energy effective 
model for iron-pnictides, one may ask how does this 
model compare to the existing models for iron pnictides 
and what are its advantages relative to other models. 

The effective itinerant models for iron-pnictides can 
be broadly divided into two groups, those which use or- 
bital or Wannier states as their basis, and those which 
use the exact eigenstates of the Hamiltonian, i.e., bands, 
as their eigenstates. In an earlier subsection, we have al- 
ready discussed the advantages our effective model has in 
comparison to the band-effective models . Among other 
issues, a model which uses the eigenstates of the Hamil- 
tonian as its basis suffers from the non-analyticity of the 
basis at the M-point. This may introduce problems in 
dealing with the states in the vicinity of that point, as 
well as issues with the interaction terms. 

The models which use orbital or Wannier states dif- 
fer by the number of orbitals/states per iron atom in the 
model. The simplest ones are two- and three-orbital mod- 
els and we discuss these in the following subsubsections. 
Speaking of models with four, five, or more orbitals, one 
can say that they reproduce the band structure of iron- 
pnictides with a great accuracy and preserving the proper 
symmetries. In fact, we rely on five orbital inodels for es- 
tablishing the symmetry properties of the states at the 
high symmetry points and also as a source for the numer- 
ical values of our model. The problem with these models, 
however, is that they carry too much information which is 
often irrelevant, but adds to the complexity of the prob- 
lem. Particularly in weak-coupling problems, one should 
not have to worry about degrees of freedom which are 
far above or below from the Fermi level. Within tight- 
binding models, removing these degrees of freedom, as 
we shall see, cannot be done without an impact on the 
qualitative features of the spectrum. 



1. Two-orbital models 

The first and the simplest effective model for iron- 
pnictides, was introduced in Ref. [S^- This two orbital 
model has been preferred due to its appealing simplicity. 

The construction of the model is based on the obser- 
vation that both hole and electron states at the Fermi 
surfaces of iron-pnictides have the highest overlap with 
the dxz and dyz states. The model uses only these two 
orbitals as the basis, and with only the nearest and next- 
nearest hopping processes between iron atoms, it is able 
to reproduce the characteristic Fermi surfaces of iron 
based superconductors: two hole pockets around the F- 
point and two electron pockets around the M-point, once 
the folding of the Brillouin zone, due to the doubling of 
the proper unit cell, is taken into account. 

The entire band structure of this minimal model for 
iron based superconductors is given in Ref. l59l . Using 
the parameters given in that paper, we reconstructed the 
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result and in Fig. [9] plotted the band structure along the 
E = r — M line for comparison with the five band mod- 
els, Fig. m The two bands crossing the Fermi level near 
the F-point correspond to the hole bands, the two bands 
that cross the Fermi level near the M-point correspond 
to the electron bands. Clearly, this model can reproduce 
the symmetry of neither hole nor electron bands. For 
example, while one of the hole bands originates from the 
iJg-doublet at the F-point, and has A2 symmetry prop- 
erties along this line, the other state originates from the 
wrong, Eu, doublet, and accordingly has Bi and not B2 
symmetry, as it should. Since Bi states are even under 
the {tr^l^^} mirror, any states on this Fermi surface are 
even under this mirror, i.e., belong to A' irreducible rep- 
resentation of Pk- Similarly, both electron bands near 
the M-point originate in _EM2-doublct. While one band, 
B2, has the proper symmetry along the S-linc, the other 
band transforms as Ai, not Bi as it should. 

More recenth;, another two-orbital model was pre- 
sented in Ref. |60- This model uses the same states as 
Ref. 153 . but it includes hoppings up to third nearest iron 
neighbor. We believe that this model violates mirror re- 
flection symmetry about the xz- and yz- planes. Such 
reflection symmetries force ii^ = —tiy and t^x = t:^,, , im- 
plying tis = and tsd = in Eq. (5) of the Ref. [60; such 
symmetry has not been reinforced as can be seen from 
the explicit values stated on page 6 of Ref. mX 
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FIG. 9: The band structure of the two-orbital effective model 
for iron-pnictides along the E = F — M line based on Ref. 
I59I . Comparison with the Fig. |4] shows that the model does 
not reproduce the symmetry for all bands crossing the Fermi 
level. 

Can any modifications of the two-orbital model recover 
the correct symmetry for its bands? We would argue that 
the answer is no — this model does not have sufficient 
ingredients to match all the bands parities correctly. The 
first problem is that it does not include the ii^Ms-doublet 
states at the M-point, inevitably missing states which in 
more accurate band structure calculations describe the 
electron bands. Further, the parity of the bands under 
the {c^l, 5} mirror does match: there are two Fermi sur- 



faces with odd, and two Fermi surfaces with even states 
under the mirror. A more accurate model must repro- 
duce one electron Fermi surface which is even under the 
{cr^l^i} mirror, and the three remaining Fermi surfaces 
(two hole and one electron) which are odd under this mir- 
ror. We therefore conclude that any two-orbital model is 
unable to reproduce the symmetry properties of Fermi 
surface states in iron-pnictides. 



2. Three-orbital models 

The next in complexity are the three-orbital models. 
The two representative three-orbital models we discuss 
here were introduced in Refs. lo and l62l. In both of these 
models an additional iron orbital, dx2-y2 is included in 
the basis of the effective model. This is motivated by 
the ab initio calculations which demonstrate that the 
electron Fermi surface states have high overlap with the 
Bloch states composed of the da;2_y2 orbital. 

In Fig. [TO] we plot the dispersions for these two mod- 
els along the E-line. We notice that both models are 
able to properly reproduce the symmetry properties of 
the bands crossing the Fermi level: the two hole bands 
originate in the _Eg-doublet at the F-point; the electron 
bands near the M-point originate from -Emi- and i?M3- 
doublets. Nevertheless, each of these two models is un- 
able to reproduce the qualitative aspects of the iron- 
pnictides band structure. 

The band structure of the model from Ref. [6l| con- 
tains one additional band which crosses the Fermi level, 
thus leading to one excess Fermi surface. The authors 
acknowledge this, and conclude that, in order to remove 
this spurious Fermi surface, the model must include at 
least one additional orbital. We find these statements 
to be in agreement with the symmetry properties of the 
bands. Since the band which creates the spurious Fermi 
surface transforms according to Ai along the S-line, it is 
sufficient to include either dxy orbital in the model. One 
of the bands corresponding to this orbital transforms ac- 
cording to the Ai irreducible representation along the S- 
linc. Since the band which produces the spurious Fermi 
surface also transforms according to the Ai, the two 
bands must hybridize. At the F- and M-points, where 
the hybridization is absent, the dxy state has a nega- 
tive and positive eigenvalue, respectively. Therefore, for 
large enough hybridization the spurious Fermi surface is 
removed, but the price paid is the introduction of a new 
orbital. 

By cleverly arranging the hopping amplitudes in the 
three-orbital model, Ref. |63 arrives at the band structure 
where the spurious Fermi level is absent. By looking at 
the dispersion in Fig. [TO] we find that this is achieved 
by moving the energy of the i?„-states at the F-point 
below the Fermi level, while changing the energy of a 
single Big state from negative to positive. In addition, 
for the parameters given in Ref. l62l . this model makes the 
energy of the £'iv[2-doublet higher than that of the i?M3- 
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doublet. Thus, this model violates the energy ordering 
of the states at both F- and M-points. One can easily 
conclude, by looking at Fig. [101 that any three-orbital 
model must either a) have a spurious Fermi surface, or 
b) have the states at the F-point ordered incorrectly. 
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FIG. 10: The band structure in three-orbital eflFective models 
for iron-pnictides along the E = F — M line, reconstructed 
from Refs. a) |6l| and b) |62. The symmetry properties of the 
hole and electron Fermi surfaces states are correctly repro- 
duced in both models, however, in a) a spurious Fermi surface 
is present; in b) there are no spurious Fermi surfaces, but the 
ordering of the bands at the F- and M-points is incorrect. 



By the number of degrees of freedom, i.e., states per 
iron atom, the three-orbital models are equivalent to our 
low-energy effective theory. Any model which would in- 
clude additional orbitals in order to remedy the short- 
comings of the three-orbital models is, by construction, 
more complicated than the model we introduced. There- 
fore, we do not need to compare such models to the 
low-energy effective model we introduce, and can safely 
assume that our model is the simplest model for iron- 
pnictides which respects the symmetries of the states 
near the Fermi level. 



C. Electron-electron interactions without the 
spin-orbit interaction 

In this part of the section we present the term in 
the Hamiltonian which describes the contact electron- 
electron interaction. Such interaction may lead to the 
symmetry breaking orders, such as spin-density waves 
or a superconducting pairing. How the interaction terms 
lead to the broken symmetries is beyond the scope of this 
paper and is left for the future. Nevertheless, for com- 
pleteness, we present the corresponding terms given that 
they follow naturally from the symmetry classification. 
We focus only on the quartic term in this paper, but 
following similar steps, one can have a straightforward 
derivation of higher order terms if these are desired. 

Focusing on the quartic term, the translational invari- 
ance implies that it can be split into two parts. 



"Hi, 



H., 
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H: 
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(28) 



r(0) 



The first term, -ffj„j , contains terms that are products 
of two charge neutral bilincars, each carrying no mo- 
mentum, therefore representing small momentum trans- 
fer processes. In the second term, H^^^ , each term is a 
product of two bilinears, both having momentum M, and 
it describes interaction processes where large momentum 
M is exchanged. The Umklapp processes are contained 

m iii„t . 

A product of two bilinears has the symmetry properties 
given by the product of the individual representations. 
The product table for the irreducible representations of 
Ti^h implies that our quartic terms — a product of two 
bilinears — transforms under the trivial representation 
only if both bilinears belong to the same irreducible rep- 
resentation. The time reversal invariance similarly im- 
plies that bilincars must be either both even or both odd 
in order for the product to be time reversal invariant. 

Given a spin singlet bilinear X^o- V'J(r)F^'" i/'crli")! the 
transformation of V''s under the space group operations 
and the time-reversal induce the transformations on 
VYj ■ The 36-dimensional representation, given by the 
F's has been reduced and presented in Appendix [Cl For 
one-dimensional irreducible representations m ^ \. For 
the two-dimensional ones, to = 1 or m = 2 denoting 
the component of the irreducible representation. Two 
irreducible representations, A\^ and i?^, correspond- 
ing to z = 1 and i = 7 respectively, have multiplicity 
All = AI7 = 3. In other words, there are three linearly 
independent bilinears, enumerated by j = 1,2,3, that 
transform according to Af and S^ respectively. For all 
the other irreducible representations, the multiplicity is 

Mi = 1. 

r(0) 



We divide H^ J into a direct and mixed part. 
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int. mix" 



(29) 



The direct part contains products of two identical spin 
singlet bilincars. All the terms that are products of two 
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different bilinears belonging to a same irreducible rep- 
resentation are contained in the mixed part. We can 
therefore write 



12 Mi 



^ A 1 A 1 



.gp^x 



i=l J = l 
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E I E V'i(r)rgV.(r) 

m=l \o-=-f-,4. 



(30) 



for the direct, and 

r(0) 



H. 



J/dr EE E ^^ 



;;(«■')- 



i=i,4i=i j'=i+i 



(31) 

for the mixed part. 

Following the same argument, the quartic contact in- 
teraction term, iJ^^^ , reads 



H.. 



(M) 



^/dr E.9.E f E ^iWrl:V.(r) I 



(32) 

The multiplicity of each irreducible representation E-^^ 
is exactly one, hence we do not sum over j in this expres- 
sion. So far, we have identified 30 independent coupling 
constants, 16 in Eq. (gni), 6 in Eq. (|3I]), and 8 in Eq. ([5^ . 
In addition to the singlct-singlct quartic terms we have 
just presented, there are symmetry allowed triplet-triplet 
terms. For each term in Eqs. (p0| -(|32 |) . there is an anal- 
ogous triplet-triplet term, schematically obtained by the 
substitution 

where each F stands for one of the F,-*™ matrices. The 
summation is assumed only over the spin indices writ- 
ten in Greek letters. The terms obtained in this manner 
obey the lattice, time-reversal and spin SU(2) symme- 
tries. Seemingly, Eq. ([55]) introduces additional 30 in- 
dependent coupling constants, but this is not true for 
the contact interaction. The Pauli matrix completeness 
relation, 



app' = 2Sal3'Sa'l3 — Saa'Sjl^' 



(34) 



can be used in each triplet-triplet term. Rearranging 
fermion operators in such an expression yields a quar- 
tic interaction term that is a product of two spin singlet 



bilinears, Eqs. (|30|) - (|32l) . Hence, none of the 30 triplet 
coupling constants, g^*\ is independent. We conclude 
that the terms in Eqs. ([301)- (EH represent the most gen- 
eral quartic contact interaction terms, with spin SU(2) 
symmetry, in our low-energy effective model. 



IV. SPIN-DENSITY WAVE ORDERS AND 
SOME PROPERTIES OF THE BAND 
STRUCTURE IN THEIR PRESENCE 

The only spin-density wave order parameter at the 
wavcvector M with sublatticc magnetization peaked on 
the iron sites corresponds to the irreducible representa- 
tion Em4 (see Figs. [2] and [3]). The order parameters, cor- 
responding to i^Mi, £'M2, and Em.3, vanish on the iron 
sites. Since they correspond to a different irreducible rep- 
resentation, they also must have different transition tem- 
peratures unless, of course, the system is fine tuned to 
a multicritical point. That's because there is no symme- 
try which guarantees the coefficients of the second order 
invariants to vanish at the same temperature. 

When the spin SU(2) symmetry is present, the irre- 
ducible representations discussed in the Section |ll] can 
be used to classify the orbital component of the spin- 
density wave order parameters; the spin magnetization 
component is independent. If the spin-orbit coupling is 
included the spin SU(2) symmetry is absent. In such a 
case, under a symmetry operation, the direction of the 
spin must be transformed together with the spatial coor- 
dinates. The spin-density order parameters can be clas- 
sified according to the irreducible representations of the 
space group determined in Section |ll] as long as we are 
careful to include the sublattice magnetization when per- 
forming the symmetry operations. Any spin-density wave 
order parameter is odd under time-reversal. 

With spin-orbit coupling, the spin part of the order 
parameter transforms under the lattice space group op- 
erations together with the orbital part. Being an axial 
vector, the spin part of an order parameter transforms 
under Eqs. (|22I -[M )) . Therefore, when writing down the 
invariants, we must first determine the overall symme- 
try properties of a spin-density order parameter, A|;a . 
Here, b is the spin component index and E^^ corresponds 
to the irreducible representation of the orbital part. This 
is done by multiplying the corresponding representation, 
and the result is presented in Table 1x1 

To illustrate the connection between the orbital part 
and the overall symmetry of a spin-density wave or- 
der parameter, let us consider the collinear spin-density 
wave^^iS shown in the upper left corner of FigfTTl The 
orbital part of this order transforms as one of the i?M4 
components, let us set it to E^^. Then, the local mo- 
ments on irons are pointing in the X-direction too^i^, 
yielding the overall order parameter symmetry E^^, ac- 
cording to Table |Xl Had the spin pointed in Y- or z- 
direction, the overall order parameter symmetry would 
have been E^2 or ^m3i respectively. 
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FIG. 11: The magnetic moments on iron (black) and pnictide 
(red) atoms when only a single component of the overall order 
parameter condenses. The moments on both iron and pnictide 
sites are locked into X-, Y-, or z-direction. 
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FIG. 12: The magnetic moments on iron (black) and pnic- 
tide (red) atoms when both order parameter components for 
a particular irreducible representation condense at the same 
time. The local moments are locked into x-, y-, or z-direction. 
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Table |X] may also be used in determining the symme- 
try properties of the orbital part of a spin-density wave, 
once the symmetry properties of the overall order param- 
eter have been set. In the physically relevant case, where 
the overall order parameter is E^^, we find that, next to 
Apx order parameter, there are two other spin-density 

wave order parameters A^^ and A^ ^ which belong to 
the same irreducible representation. These order param- 
eters are automatically induced by the presence of A^x ■ 
The latter of the two has the orbital part transforming 
according to £'m2, and as such, it is finite on pnictogen 
atoms. This leads us to the prediction that the collinear 
spin-density wave of the kind reported in Refs. l32ll33l must 
induce moments on pnictogen atoms in the presence of 
spin-orbit coupling. 

For the overall order parameter £'^2' where the mo- 
ment on iron sites points in the F-direction (see Fig. [TT|) , 
there are no corresponding spin-density wave order pa- 
rameters with an £'m2 orbital part. Therefore, such a 
collinear spin-density wave induces no moment on pnic- 
togen atoms. For the overall order parameter Ej^^, where 
the moment on iron sites points out of plane, the induced 
spin-density wave results in the moment on pnictogen 
atoms pointing in-plane. In both cases where a collinear 
density wave with a moment on pnictogen atoms is in- 
duced, the orbital part of the induced order parameter is 

-^M2- 

There is one additional spin-density wave order param- 
eter, namely the one with the overall symmetry proper- 
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TABLE X: This table gives the overall symmetry properties 
of a spin-density parameter A^a , whose orbital part trans- 
forms as a-component of the E-Mi irreducible representation 
(left), while the spin points in the b direction (top). The 
table can also be used to determine the orbital part a spin- 
density wave order with an overall -Bmj symmetry when the 
spin points in a particular direction. Omitted from the table 
is the parity under the time reversal: the overall order param- 
eter is odd under the time reversal thanks to its spin part; the 
orbital part is even. 



ties given by E-M4- A spin-density wave order parameter 
with such an overall symmetry cannot have an E-M.i or- 
bital part, as seen in Table [Xj therefore implying that 
such a spin-density wave-order parameter carries no mo- 
ment on the iron sites (see Fig fTTj) . 

The spin-density wave orders which we consider is ei- 
ther a collinear spin-density wave or a co-planar four-fold 
symmetric density wave shown in the upper left corners 
of Figs. fTTj and \T% The collinear state, with the sublat- 
tice magnetization peaked on iron, occurs when only a 
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single component of an -Emi order parameter condenses. 
This state is ubiquitous in iron pnictides^i^^. The four- 
fold density wave may occur when both components of 
an E-Mi order parameter condense with the same magni- 
tude. Such a state may have been observed recently^. 

From the symmetry properties of each component of 
the low-energy effective theory spinor in Eq. (|17|) . we 
know how to construct symmetry breaking terms. Since 
the magnetic moment of these orders is situated on iron 
sites, the orbital part in each symmetry breaking term 
is a bilinear which transforms as one of the -EM4 compo- 
nents. These are constructed and analyzed in the follow- 
ing subsections. But first, we address the issue of double 
degeneracy of all bands, as found in the dispersions with 
the spin-orbit interaction absent. Fig. |6l or present. Fig. 
[51 in the model. 



A. Kramers degeneracy 

In a model with no spin-orbit interaction, each elec- 
tron state is doubly, or even-, degenerate. This is 
an automatic consequence of the spin SU(2) symme- 
try: the Hamiltonian is diagonal in the spin space, 
hence for each Bloch state with spin up, V'kt(r) = 
1 



„2kr 



"k(r) 







there exists a Bloch state with the 



same orbital part, Uk(r), but with the opposite spin, 




V-kilr) 



■"uklr) 



The two states are orthog- 



onal, and, since they arc related by a symmetry of the 
Hamiltonian, they have the same eigenvalue Ck- 

Once the spin-orbit interaction is present in the prob- 
lem, the SU(2) symmetry is lost and the double degener- 
acy of each state does not follow automatically. It may 
or may not exist, depending on the remaining symmetry. 
Generally, if the system is left invariant under an inver- 
sion followed by the time-reversal, for each single electron 
eigenstate of the Hamiltonian at k, there is another one 
which is orthogonal to, and degenerate with it, at the 
same k. Note that, for such Kramers pairs to exist, it 
is required that only the product of the inversion and 
the time-reversal is a symmetry of the system; the two 
operations individually need not be good symmetries. 

To show that this holds, even when the inversion is fol- 
lowed by a translation, as is the case here, let us assume 
that one eigenstate of the Hamiltonian at k is 



*k(r) 



Wk(r) 
Wk(r) 



(35) 



Since the Hamiltonian is, in general, not spin SU(2) sym- 
metric, an eigenstate will have both of its spin-up and 
spin-down components non-zero. The Kramers partner 
of ^k(r), obtained by the product of inversion and time- 
reversal (assumed to be a symmetry of the Hamiltonian), 



vl/k(r) = e{^\T}^u{v) - ee"^(---) ( ^^^ ^J 



^k-r^-ik-r / Vy,(T - r) 



-zk-T 



(icr2)A'*k(-r-r). 



(36) 



The inversion operation followed by a translation by vec- 
tor X, {^It"}, is equivalent to an inversion with respect 
to the point at t/2. Any inversion keeps the spin part 
of the wave-function invariant, as it must, because spin 
is an axial vector. The time-reversal, Q, complex conju- 
gates the value of the wave- function (K) , while acting as 
i(j2 in the spin space. Since ^k(r) is a Bloch state with 
momentum k, from Eq. (p6)) it follows that ^k(r) is also 
a Bloch state with the same momentum. 

Under the assumption that 0{i|T} is a symmetry of 
the Hamiltonian, it follows that the Kramers partner Eq. 
([M]) of the state ^k(r) has the same eigen-energy. These 
two states are orthogonal. 



(37) 



dr *;;(r)*k(r) = 0, 



which follows directly from the definitions in Eqs. ()35|) 
and ([55]). 

In the absence of a symmetry breaking order, any op- 
eration {i|t-|-To}, where t is an integer translation vector, 
is the symmetry of the Hamiltonian. Unless there is an 
externally applied magnetic field, the time-reversal sym- 
metry is also automatically present, with or without the 
spin-orbit interaction. Therefore, the product of these 
two operations is also a symmetry, implying the Kramers 
degeneracy. This is in agreement with the dispersions 
shown in Fig. [51 



B. CoUinear spin-density wave order parameter in 
the low-energy effective model 

Since the magnetic moment rests on iron atoms, the 
orbital part of the coUinear spin-density wave order pa- 
rameter transforms according to one of the components 
of the E-M.4. Accordingly, this order parameter couples to 
any bilinears in the low-energy effective theory for which 
the orbital part has the same symmetry properties. The 
only two bilinears in the low-energy effective theory form- 
ing an orbital £'m4 doublet are 



E^4-- E^r,Jk) 

k 

^M4: E^r,a(k) 



1 



-1 



CTa/3V'y,/3(k) + /i.c., (38) 
o^Q^V'x,/?(k) +h.c..{39) 



One of these two bilinears appears in the symmetry 
breaking term in the presence of the spin-density wave 
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order. Wc chose the colhnear spin-density which is anti- 
periodic in the X-direction, thus the orbital part of the 
order parameter is Eq. (|38| . For such a choice of the 
orbital part of the order parameter, the experiments 
find the magnetization to point in the X-direction22i^. 
The symmetry breaking term corresponding to such a 
collinear spin-density wave is 



H< 



SDW 



k Q,;3=t,i 



1 




'^i/i^rA^) + h-c 



(40) 



In Eq. (|40p . we have written only the k- independent 
symmetry breaking term of a collinear spin-density wave 
state. Wc neglect all k-depcndent terms because the radii 
of the hole and electron pockets are relatively small com- 
pared to the extent of the Brillouin zone. 

In Fig. [T^ we plot the dispersion of the low-energy 
effective model with the symmetry breaking term, Eq. 
(PO)) . The plot shows both the result in the absence of 
any spin-orbit interaction (dashed lines) and with the 
spin-orbit coupling terms, Eqs. (P5|) and ([?7)) . included 
(solid line) . The value of the spin-orbit coupling is set to 
Ar = Am = 80meV according to Ref.[2^ The value of the 
symmetry breaking parameter is set to Asdw = 65meV 
since this value produces Fermi surfaces similar in shape 
to those found experimentally^iH. 

All bands in the dispersion in Fig. [T2] are doubly de- 
generate, and find this to be true even away from the 
high symmetry lines. This is the Kramers degeneracy, 
described in the previous subsection. The collinear spin- 
density wave state is odd under the time reversal; it is 
also odd under any inversion {i\t + W}, where t is an 
even translation. Therefore, the collinear spin-density 
wave is even under the product of these two operations 
and all the states come in Kramers pairs. 

In the absence of any spin-orbit coupling, the disper- 
sion of one of the electron bands is not affected by a pres- 
ence of the collinear spin-density wave. In our example, 
it is the two ipx.u components of the spinor, Eq. ()17p . 
that have are decoupled from the remaining degrees of 
freedom, even with the spin-density wave term, Eq. (|40| , 
present. Had we considered the collinear state where the 
other component of the i?M-doublet condenses, it would 
be the two components in ij}Y,a that are decoupled and 
unaffected by the presence of the spin-density wave order. 
Either way, there arc two bands decoupled from the oth- 
ers and the crossings of these two bands with the others 
results in line degeneracies. These crossings are guaran- 
teed by the {cr^jj^} nurror reflection symmetry. Even 
if one considered a spin-density wave order parameter 
which is kx- and fcy-dependcnt, these lines would remain. 
However, if the spin-density wave order parameter has a 
/cz-dependence, such that it has a piece which is a k^ odd 
function, then, based on the symmetry, additional terms 
would appear in Eq. (|40[) which would lead to avoided 
crossings between the bands. When the moments are 




b) 



FIG. 13: The band structure from the low-energy effective 
model in the presence of a collinear spin-density wave or- 
der, a) The dispersion of the low-energy effective model 
along the high symmetry lines. The orbital part of the 
spin-density wave order parameter transforms as the -Em4- 
component. The spin points in the X-direction. Without 
the spin-orbit coupling (thin dashed lines), both Dirac nodes 
and lines where bands cross are found in the dispersion. The 
spin-density wave order splits the four-fold degeneracy at the 
M-point. With the spin-orbit coupling, all the band crossings 
are avoided, although the splittings are too small to be visible 
in the plot. Each state is doubly degenerate, a consequence 
of the Kramers doublet existence in the collinear spin-density 
wave phase, b) The Fermi surfaces bear resemblance to the 
experimentally determined Fermi surfaces^i^. 



localized on iron sites only, as it is in the absence of spin- 
orbit coupling, the spin-density wave order is always an 
even kz function, hence the Dirac lines are unavoidable. 
This can also be shown following the Wigner-von Neu- 
mann type of arguments for the presence of degeneracies 
in the system. 

In our example in Fig. I13[ the Fermi surface crossing 
band which is decoupled from the hole bands is the inner 
electron band in the right hand side of the plot {Bi in 
Fig. SJi, and Fig. [5|) , and the outer electron band (also Bi 
on Fig. [5]) on the left hand side of the plot. For a k not 
on a high symmetry line, the electron band not affected 
by the spin-density wave term Eq. (|40p is the A' band 
in Fig. [5l On the other hand, the electron band coupled 
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by Eq. (gO]) to the hole bands is the A" band in Fig. [H 
which also has B2 symmetry properties along the high 
symmetry lines plotted in Fig. [131 

Putting the two decoupled "spinor" components, ipx.cr, 
aside, we note that the collinear spin-density wave term 
in Eq. (|40)) couples the first components of ipr,a, and the 
first component of 'ipY,tT- This leads to a gap opening be- 
tween the hole and electron bands for k's which do not 
lie on the high symmetry lines. For k's on the high sym- 
metry lines, plotted in Fig. I13i we notice that only one 
band crossing, on the left hand side of the plot {kx — 0) 
between the lower hole band and one electron band, is 
avoided; all the other band crossings remain in the ab- 
sence of spin-orbit interaction regardless of the strength 
of the spin-density wave order parameter Asdw hi Eq. 
(Uni) • Had we considered a k-dependent spin-density wave 
order parameter AsDw(k), a gap would open at another 
band crossing, the one on the right hand side of the plot 
{ky = 0) between the lower hole band and one electron 
band. 

Apparently, the spin-density wave order does not open 
a gap between the upper hole band and the electron band 
along these two directions, thus leading to Dirac cones 
in the dispersion of the spin-density wave state in iron- 
pnictides^. These crossings arc protected by the vertical 
mirror reflection {cr^lH} on the ky — line (and simi- 
larly by {(T''^ I i I } on the kx = line) . When we classified 
the bands along the S-line, the two hole bands, origi- 
nating from the £^g-doublet at the F-point, transformed 
according to A2 and B2 irreducible representations, re- 
spectively. Therefore, one hole band is even, and one 
is odd, under the mirror reflection {c^|^^}. The sym- 
metry breaking term Eq. (j40p couples the electron B2 
electron band along this line to only one of the two hole 
bands. In our example, it is the B2 hole band. If there is 
any coupling between the electron band and the A2 hole 
band, it must be odd in ky so that the overall symmetry 
breaking term has the same parity under the {c^l^^} 
mirror reflection as the ky even term which couples the 
B2 hole band and the electron band. The presence of 
these nodes can alternatively be deduced from Wigner- 
von Neumann type of arguments as we do in Subsection 
CVDl 

These Dirac nodes arc protected only as long as no 
other bands enter the effective model, i.e., as long as the 
Asdw energy scale is lower than the energy of the bands 
excluded in the low-energy effective theory. We studied 
the evolution of the full tight-binding modelii disper- 
sion in the presence of the spin-density wave symmetry 
breaking term on each iron orbital, and determined that 
the nodes vanish once the Aiu band at the F-point gets 
close to the Fermi level. This happens for the values of 
Asdw ~ lOOmeV. 

Alternatively, we found the Dirac nodes to be instanta- 
neously removed once the spin-orbit coupling is included 
in the Hamiltonian. This is a consequence of the coupling 
of the two hole bands {A2 and B2 along the S-line) in Eq. 
(j26p . Since, in the presence of the spin-orbit interaction, 



each hole band on the E-line is an admixture of both Eg 
doublet states, in a collinear spin-density wave state a 
gap opens at each crossing of the hole bands with the B2 
electron band. Similarly, the spin-orbit coupling also cou- 
ples the two electron bands, thus it leads to the removal of 
the all the degeneracy lines where the dispersion of tpx.a- 
states used to intersect with the other dispersions in the 
absence of the spin-orbit interaction. Therefore, we flnd 
that the dispersion of a collinear spin-density wave state 
has no four-fold degeneracies anywhere once the spin- 
orbit coupling is present in the model. The four-fold de- 
generacy at the M-point, which was protected even with 
the spin-orbit coupling Eq. (p7)) . is also split due to the 
collinear spin-density wave order. 



C. The four- fold symmetric (coplanar) spin-density 
wave in the low-energy effective model 

In a four-fold symmetric coplanar state, in addition to 
one component of the order parameter, here £"^1' ^^^ 
other component of the total order parameter, E^^, also 
acquires an expectation value, which is the same in mag- 
nitude. The symmetry breaking term, therefore, con- 
tains a term equivalent to Eq. (|40p . and its E-m doublet 
partner, i.e., the term obtained by the action of {cr^jOO} 
mirror reflection. 






4.J^) 



1 




4J^) 




-1 



'^Ipi'xA^) 



+ h.c. (41) 



In Fig. [13] we plot the dispersion in the low-energy ef- 
fective model with this symmetry breaking term present. 
We choose the strength of the order parameter to be the 
same as in the previous subsection, Ac4 = 65meV. The 
first notable property of the dispersion is that the bands 
are not doubly degenerate, i.e., the Kramers degeneracy 
is broken. One can understand this by noticing that there 
are two symmetry breaking terms in Eq. (j4ip . For any 
lattice inversion, {i|t -I- Tq}, where t is an even lattice 
translation vector, the first term in Eq. (|^T|) is odd un- 
der the action of this lattice symmetry; the second term 
is even. If we consider the inversions where t is an odd 
lattice translation vector, then the first term is even and 
the second term is odd. Since both terms are odd under 
the time reversal symmetry operation, the entire sym- 
metry breaking term is neither even nor odd under the 
combined transformation, therefore, for an arbitrary k 
the Kramers degeneracy is lifted as seen in Fig. [TD 

Since the Kramers degeneracy is broken in the four- 
fold coplanar spin-density wave, for an arbitrary k we 
find each band to be non-degenerate. However, at F and 
M we find each state to be doubly degenerate. The de- 
generacy is exact only at these points, while for arbitrary 
small momenta, the dispersion is linear, i.e., we find six 
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FIG. 14: The band structure of the low-energy effective model 
in the presence of the four-fold symmetric spin-density wave 
order. The local moments are plotted in Fig. [121 ^) The 
dispersion of the low-energy effective model along the high 
symmetry lines without the spin-orbit coupling (thin dashed 
lines), and with the spin-orbit coupling (solid lines). The 
Kramers degeneracy is broken at every k-point. On the right 
hand side of the plot {kx = or fcy =0) there are no band 
crossings. On the left hand side of the plot (k^ = or fcy = 0) , 
the bands crossing is allowed leading to Dirac points in the 
spectrum of this phase, b) Thanks to the absence of the 
Kramers degeneracy, the Fermi surfaces appear in weakly split 
pairs. 



Dirac cones. Wc also find other Dirac points in the spec- 
trum along the two special directions, k^ = or ky = 0. 



D. Wigner-von Neummann analysis of the 
degeneracies: Dirac points and lines 

Until now, we have analyzed the excitation spectrum in 
the non-superconducting states of iron based supercon- 
ductors by diagonalizing particular model Hamiltonians. 
Such an approach is useful as it provides us with an in- 
tuition about the presence or absence of degeneracies in 
the spectrum. More generally, the degeneracies, even the 
accidental ones, are deeply tied to the symmetry of the 
problem. In this section we use a variant of the Wigner- 
von Neumann argument to prove when the degeneracies. 



such as Dirac nodes or lines, may appear in the spectrum, 
and when the bands must avoid each other. 

The basic idea is to consider two variational states 
I^'Ij. ) and 1^]^. ) for each k. Each of these states is 
assumed to be symmetry adapted and chosen to trans- 
form according to one of the irreducible representations 
of the group of the wave- vector at k, i.e., Pk. Then, 
the Hamiltonian is projected onto these states for each 
k. The symmetry properties of these states dictate the 
form of the projected Hamiltonian which can, in turn, be 
used to determine the feasibility of the band touching, or 
whether additional fine tuning is necessary. 

We are mainly interested in the "high-temperature" 
state, which breaks no symmetries, and the collinear 
spin-density wave state. For both of these states, each 
band is doubly degenerate due to the Kramers theorem, 
as proven in the previous subsection. Therefore, an as- 
sumption that the two bands cross requires us to consider 



at least four states: |^_ 
partners \^i^ ), and \Wj^ 



(1) 



I* 



(2) 



and their Kramers 



The coordinate and spin rep- 
resentation for these states are given in Eqs. ([55)1 and 
36|) . In the basis formed by the four states. 



(i<^^),i^L^\i4^^),i*L'') 



(42) 



the projected Hamiltonian is a 4 x 4 matrix 

3 

Hf = a{k)T3 (g)l + 6(k)ri ® 1 + ^ c,(k)r2 ® a^. (43) 

Seemingly, there are sixteen independent coefficients, but 
the fact that I'I'k ) and 1^^ ) come with their Kramers 
partners, reduces the number to five. This can be seen by 
observing that under the inversion and the time reversal. 



e{i\T}Hf = (1 ® fT2) [KHf] (1 ® CJ2) 



(44) 



When the spin-orbit interaction is absent, the Hamil- 
tonian has the spin SU(2) symmetry, and any spin rota- 
tion is independent of the lattice symmetry operations. 



,(1 



,(2) 



Both l^j^ ') and \^\^ ) can be chosen to be spin up, and 
their Kramers partners must be spin down. With such a 
choice, all three coefficients Ci(k) must vanish for every 
k. One is then left with only two independent coeffi- 
cients, a(k) and 5(k), in Eq. (|^5)) . Therefore, in order 
for the degeneracy to appear, it is enough for the two 
functions to vanish at the same k. Since, for k^ = 0, 
each is a function of two variables, they vanish along one- 
dimensional curves. Two such curves may generically in- 
tersect without any further fine tuning. This would lead 
to Dirac points, provided that |$j^ ') and j^pj^ ') have the 
same parity under the in-plane mirror rcfiection followed 
by a 'half-translation'. If the two states have an oppo- 
site parity, then, the Hamiltonian must commute with 
T3 (8) 1, and this sets 6(k) = for all k's. We see that 
in such a case, only a single function of two variables 
a{kx, ky) must vanish in order for accidental degeneracies 
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to appear. Since this may happen along one-dimensional 
curves without any further fine tuning, the two bands 
may intersect along lines in k-space. In fact, we may al- 
ready see this in Eq. ([TB)) . where /i^(k) and h^Ck) are 
not coupled; this holds to any order in k because, one 
block corresponds to the states even and the other block 
to the states odd under the in-plane mirror. This, in prin- 
ciple, allows one to 'unfold' the two-iron Brillouin zone to 
one-iron Brillouin zone; as soon as spin-orbit is taken into 
account, such unfolding can never be done. The above 
argument goes through in the collinear spin-density wave 
state because one can always choose a 'half-translation' 
so that such a state is invariant under the combined op- 
eration with the in-plane mirror. Additionally, vertical 
plane mirror operations can be used to locate the posi- 
tions of the Dirac nodes along the Brillouin zone diago- 
nal, thus making connection with the results of Ref. |34| . 
However, as we will now show, including the effects of the 
spin-orbit coupling in the collinear spin-density wave re- 
moves all the accidental degeneracies without additional 
fine tuning. 



1. Absence of degeneracies at a low symmetry k 

For an arbitrary k ~ {kx,ky,0), away from the high 
symmetry points and lines, there is only a single sym- 
metry operation in P^. In the high temperature state, 
where no symmetries are broken, the mirror reflection, 
followed by a 'half-translation' {cr^||^} corresponds to 
such an operation. On the other hand, in the collinear 
spin-density wave, the unit cell is doubled. The symme- 
try breaking term in the Hamiltonian may be even or 
odd under {cr^jg^}, depending on whether the overall 
spin-density wave order parameter transforms according 
to i^Mi O'^ ^Mii respectively. If the symmetry break- 
ing term is even under {cr^l^nlj then it is odd un- 
der {c^lj^}' ^^^'^ ^^'^^ versa. Therefore, when analyz- 
ing the degeneracies in the collinear spin-density wave, 
only the one which leaves the Hamiltonian invariant, i.e., 
the 'even' one must be considered. To keep the discus- 
sion general, therefore we assume that {ct^Itz} leaves the 
Hamiltonian invariant, where Tz will be specified accord- 
ing to the state we analyze: E^^, -Emii '^^ the high tem- 
perature symmetric state. 

l^^l- ) is either even or odd under the {ct^Itz}; simulta- 
neously, the parity of its Kramers partner must be oppo- 
site. Similarly for l^j. ). Therefore, the Hamiltonian in 
Eq. (j43|) must commute with either 1(8) (72 or with T3(E)<T2, 
depending on whether the parities of \^\^ } and I^J^ ) 
are the same or the opposite. In the first case, ci(k) 
and C2(k) must vanish for every k; in the second case, 
6(k) and C3(k) must vanish. In either case this requires 
three functions of two variables to vanish simultaneously, 
which does not happen without fine tuning. 



2. Absence of degeneracies on the F — M line (excluding 
the high symmetry points) 

If we consider the effect of additional symmetries at a 
generic k-point along the high symmetry line E = F — M 
excluding the points F and M, then there are additional 
constraints on the functions in Eq. P5|) . Note, however, 
that the constraints also remove one independent vari- 
able, because, unlike in the above discussion, we must 
move along the high symmetry line. 

The corresponding symmetry operation is {cr^|Ty}. 
Here, analogous to the above, ly- is the 'half-integer' 
translation which, combined with the (vertical) mirror 
along the face diagonal leaves the Hamiltonian invariant, 
even in the collinear spin-density wave state. When act- 
ing on our basis Eq. (|^ . up to an overall phase factor, 
it can be represented either by 1 (8 (72 or by t^ (SD cr2, 
depending on the relative parities of u{^ (r) and u\^ ' (r) 
under {cr^|Ty}. If ci(k) = C2(k) ~ 0, as required by the 
in-plane mirror, then, either C3(k) = or b(k) = 0, re- 
spectively. In either case, we are left with two functions 
of one variable, which do not vanish at the same point 
without additional fine tuning. If 6(k) = C3(k) = 0, then, 
either ci(k) = or C2(k) = 0, respectively. Again, in ei- 
ther case we are left with two functions of one variable 
which do not vanish simultaneously unless fine tuned. 
Therefore, whether in the normal state or in the collinear 
spin-density wave state, the bands generically avoid each 
other. This effect is captured within our low-energy ef- 
fective Hamiltonian. 



3. The high symmetry points T and M 

Looking at the Table U we see that, at the F-point, 
if the symmetry adapted functions come from one di- 
mensional representations, then there is no constraint on 
a(k). Since, there is no longer any freedom in choos- 
ing k, it being at F, the bands avoid each other without 
additional fine tuning. 

If the symmetry adapted functions come from a two- 
dimensional irreducible representations, we should note 
that, because the two sates must have equal parity under 
the in-plane mirror, only the case ci(k) = C2(k) = 
needs to be considered. Similarly, the parities under the 
vertical mirror {(7^|tv} are opposite, therefore &(k) = 0. 
The remaining generator {(T^|00} requires the effective 
Hamiltonian to commute with ti (E) {(Ji — (72)- While this 
forces a(k) = 0, it puts no constraints on C3(k) at F. 
Again, since there is no longer any freedom in choosing 
k, the bands avoid each other without the fine tuning. 
This effect can be clearly seen in Fig. [8^. Note that the 
remaining term has a form analogous to the Kane and 
Mele term in graphene^. 

Because, as long as no symmetries are broken, at the 
M-point there are only two-dimensional representations, 
we should consider e^^''^u\J{r), and eJ^''^uj^{r) trans- 
forming as a doublet of the same irreducible represcnta- 
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tion. Note that this is not the case in the cohinear spin- 
density wave. The two components of any E-^i doublet 
have the opposite parity under the in-plane mirror, hence 
&(k) = C3(k) = 0. Under the vertical mirror {cr^lry}, 
the two components of -Emi doublet have the same par- 
ity. The same is true for the E-M2 doublet. Therefore, 
for -Emi and £'m2, ci(k) = 0. Conversely, for Ems and 
E-Mi, the two components of the doublet have the oppo- 
site parities under {cr^|Ty}, giving C2(k) =0. The verti- 
cal mirror {cr^ |00} requires the effective Hamiltonian to 
commute with ti ® (cri — CT2)- Similarly, {ct^|00} requires 
the effective Hamiltonian to commute with ti ® (ci +(72). 
Therefore, it must commute with both ti(8)CTi and Ti(8)cr2. 
This automatically sets a(k) = 0, as well as, the remain- 
ing C2(k) = or ci(k) = 0. Therefore, the two Kramers 
degenerate bands touch at the M-point. Moving away 
from it, splits it linearly in momentum as seen in Fig. 

In the presence of the coUinear spin-density wave, the 
symmetry is broken such that the unit cell doubles and 
the M-point folds to F. Our discussion on the absence 
of the degeneracies at the F is applicable; therefore, the 
two Kramers bands avoid each other in the absence of 
fine tuning. For the value of the coUinear spin-density 
wave order parameter used in Fig. [T31 the splitting of the 
lowest, Ems, doublet is ^ O.lSmeV, and for the other, 
Em.1, doublet is ^ 7.8mcV. 



V. SUPERCONDUCTIVITY 

In this section we describe the superconducting state 
using our low-energy effective theory introduced earlier. 
We will only consider pairing with no overall momen- 
tum. This is consistent with the current experimen- 
tal findingsi^iiii^; adding finite momentum pairing is 
straightforward. The pairing term therefore involves ei- 
ther two states at F or two states at M, allowing us to 
analyze them separately. 

In the first part of the section we consider only spin 
singlet pairing. All such pairing terms in the low-energy 
effective theory can be classified according to the sym- 
metry of the lattice, Pi/nmm. We show that the Aig- 
superconductivity can be well described by three param- 
eters. Depending on the relative sign and size of these 

parameters, the ground state of our model is s++-, s_| -, 

or a nodal-s-wave superconductor. This effective model 
successfully describes the gap anisotropy for the electron 
Fermi surfaces near the M-point. In the absence of spin- 
orbit coupling, the gap is isotropic on the hole Fermi 
surfaces near F. In the second part of this section, we 
show that spin-orbit coupling necessarily admixes a spin 
triplet component which can give rise to an anisotropy of 
the gap near the F-point. 

We begin by constructing spin singlet pairing bilinears: 



where Mab must be a symmetric matrix, and ipk a' ^'^'^^ 
a = 1, ... 6 are the components of spinor Eq. ([T7)) with 
spin label a. The form of the matrix Mat determines the 
symmetry properties of the pairing term, Eq. (j45|) . Since, 
for any spin projection, the individual components of the 
spinor Eq. ([TT]) transform according to Eg, i^Mi, or Ems, 
the symmetry properties of the pairing term Eq. (|45p are 
easily deduced from the multiplication tables for these 
irreducible representations. Tables IIVI and |Vl Under a 
space group operation R we have 



Dac,bdi-^) 



C(-k)(za2)„^^^(k) 



cd L k 

(46) 

Note that, because this term is a spin singlet, the spin 
structure is left invariant under R. This is not the case 
when we consider spin triplet pairing terms, later in this 
section. 

In order to find the symmetry properties of the pair- 
ing terms in i^singiet, we decompose the outer product 
D^ (R)®D(R) into the irreducible representations of the 
group of the wave vector at F- and M-points, Pr and 
Pm • The coefficients in the decomposition are unaffected 
by the transpose due to the orthogonality of the charac- 
ters, which are, of course, invariant under the transpose. 

At the F-point such a decomposition reads 



Ea®Ea 



Alg ® Big ® B2g < 



'A 



2s- 



(47) 



Because the A2g behaves as a z-componcnt of an axial 
vector, Lz, it corresponds to an odd angular momentum 
pairing. This makes it odd under the exchange, and in- 
compatible with the spin singlet; it would require Mab to 
be antisymmetric, contrary to Eq. (P5|). Therefore, the 
pairing terms at the F-point may transform only accord- 
ing to one of the first three terms in Eq. (|T7)) . The Aig 
corresponds to the s-wave symmetry; the Big and i?2g 
correspond to dx2_y2- and dxy-wave symmetry, respec- 
tively. We focus on the s-wave symmetry. The corre- 
sponding Alg pairing bilinear reads 



^VL(-k)^r,T(k). 



(48) 



The symmetry properties of the singlet pairing terms 
at the M-point can, similarly, be determined from the 
products 

Emi ^ Emi = Em3 ® Em3 = Alg © B2g © A2U © Blu, 

(49) 

EMi®EM3 = Eg©E„. (50) 



H. 



singlet — 7 , 
k 



C(-k)(ia2)a/jMa;,V'^(k) + h.c, (45) 



In the first equation, the _Bi„ can be ignored, be- 
cause it corresponds to an antisymmetric combination 



of fcrmionic operators at the M-point, and is incompat- 
ible with the spin singlet. The two independent s-wave 
pairing terms are 



E E v^r.,(-k) 

k a=X,Y 

E E <;(-k) 



k a=X.Y 



1 





1 



^a,T(k). 



(51) 

(52) 



The remaining pairing terms, i.e., the ones not trans- 
forming according to Aig, as well as those with the (tt, tt) 
pairing momentum, are spelled out in detail in Appendix 

El 



A. Singlet spectral gaps near V and M 



25 



induced Aig spin triplet admixture. The latter will be 
considered in the following subsection. 

For the low-energy states near the M-point, the spin 
singlet pairing term with the Aig symmetry is 



H, 



M,SC 



k 



E 

i=X,Y 



i^lii-^) 



Ami 






Am3 



V'a,t(k) + h.c. 
(57) 



We define the Nambu spinor at the M-point as vl>j^(k) ~ 
(V^_^(k),V'^,^(k),^5^^(-k),^^j^(-k)). The resulting 
Bogoliubov-de Gennes Hamiltonian is 



H, 



M3dG 



k 



*M(k)^ 



^M(k) 
Am 



-/iM(k) 



*M(k), 

(58) 



Following Eq. (pS)) . the only s-wave pairing term at F where 



iJr,sc = Ar^Vf4(-k)^r,t(k)+/i.c., (53) 



/iM(k) = 



_ '"M 



hU^) 







/iM(k) 



(59) 



where we choose Ar to be real. In order to determine 
the quasi-particle energy spectrum near F, we define the 
Nambu spinor \I/r(k) = (V'r^(k), V'f ^(-k)). The result- 
ing Bogoliubov-dc Gennes Hamiltonian is 



where we used /ir(— k)* = /ir(k), which follows from Eq. 
(US. 

Notice that diagonalizing /ir(k) leaves the pairing term 
invariant. Performing a unitary transformation, lA = 

which brings us to such "band" basis, we 



U 
therefore find 



^^ I /ir(k) Ar]L2 



/ 



.(1) 



ArlL2 -/ir(k) 



where WhT{\C)U = 



U 









Ar 

V 

^ 

,(2) 



.(2) 



Ar 



(1) 



-e 
Ar 



\ 

Ar 





(155) 



The e^^'"^^ are, of 



course, the energy dispersions of the two hole bands. 
From Eq. (l55t the dispersions for the quasi-particles fol- 
low. 



E^ . 



X») 



A^ 



(56) 



We conclude that the spin singlet Aig pairing leads to an 
isotropic superconducting gap, which is identical for the 
two hole bands. Any gap anisotropy on the hole Fermi 
surfaces results either from higher order momentum de- 
pendent Aig pairing terms, or from spin-orbit interaction 



To obtain Eq. (|58|) . we have used /im(— k)* = /iM(k), 
which follows from Eq. p^ : the diagonal pairing matrix 

is Am = diag(AMi, Am3, Ami, Ams)- 

The Hamiltonian in Eq. (|58p may be brought into a 
block diagonal form such that the first block describes the 
X-components of the i?Mi and -Ems doublets, and their 
pairing, and the other block describes the ^-components 
and their pairing. The two blocks are related by the mir- 
ror reflection about the yz-plane, {cr^ |00}, and therefore 
it is sufficient to consider just the first block. This reads 



H, 



M'^,BdG 



k 



^.^^^(k),V'54(-k) 



hU^) 



6 -h+,{k) 



where 



6 = 



Ami 
A 



M3 



^x,t(k) 

V44(-k) 

(60) 



(61) 



We note that the Bogolyubov-de Gennes matrix in Eq. 
([Bn]) can be written as T3 (E) /ij^(k) +ti^ 6. Therefore, it 
anticommutes with T2I, and, as a result, for every eigen- 
statc at k with energy i?k, there is another one at k with 
energy — £'k. Therefore, the solution of the fourth order 
secular equation reduces to the solution of a quadratic 
equation only. The explicit form of the solutions is pre- 
sented in Appendix El 

Before we present the full solution, let us illustrate 
the behavior of the pairing gap on the Fermi surface 
in the limit where Ami and Am3 are much smaller 
than the energy difference between the two eigenvalues of 
/ij[^(k), e{^ ' and e{^ '. Performing a unitary transfor- 
mation which diagonalizes h^ik) = u+(k)l-|-u_(k)(T^ + 
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i'-i-(k)cr^, where 
1 



w±(k) 



k2 

-(ei±e3) + — 



1 1 

mi 7713 



-(ai ±a3)kxky 
(62) 



and 7;+ (k) was given in Eq. ((20|) , we find that the matrix 
in Eq. (jBHl) becomes 



/ 



-(1+) 







Ami cos2 ^ + Am3 sin^ ^ i (A] 
V ^(Am3 - AMi)sinxk 





~(2+) 

Ami sin^ ^ + A; 



Xk 



Ami cos2 ^ 

2(^M3 



.M3 COS^ ^ 



'k 





A* sin^ ^^ 
^M3 ^"^ 2 

sinxk 



i(A^i-A^3)sinxk \ 
A^iSin2^ + A^3COs2 2|. 





42+) 
^k 



(63) 



/ 



where cos Xk 



"-(k)/V'"-(k)^ +u+(k)2, sinxk = 

77+ (k) ± 



~(l + .2+) 
•^k 



cr(l + 



7; +(k)/^7t_(k)2+ 7;+(k)2, and 

^7i_(k)2+7,+ (k)2. 

If the band with the dispersion e)^^' crosses the Fermi 

level, while e\^ remains well separated from it, then, to 
leading order, we can neglect everything but the paring 
term which couples ej^ with —e\^ , which become de- 
generate at the Fermi level. Therefore, the gap on the 
Fermi surface — described by the zeros of ej^ 
be approximated by 

A (1 + ) A 2 Xk , A .2 Xk , 

A^ ^> =Ami cos — + Am3 sm ^ + • ■ • 



2(1+) 



can 



;(A 



Ml 



,1, A \ "-(1^) 

Am3) + ^(Ami — Am3)7(y:p) 



where the 



O 



(A: 



Ml 



corresponds to 

Vie 



?(l+)l 



"+(k) 
(64) 

terms of order 

^M3j /ICk ~ *^k I ' dropped in the 

second line. Similarly, the gap on the other band ej^ 
can be directly read off from Eq. ([M]) . however, unless 
the system is heavily hole doped, this band is always far 
below the Fermi level and need not be considered. 

Note that, unless Ami = Am3, the pairing gap at the 
M-point is anisotropic. This follows directly from Eq. 
(|63p , and can also be seen in the approximate expression 
Eq. ([M)) . Along one of the zone diagonals where k = 
(k,—k), w+(k) vanishes by symmetry. This implies that 

^k ~ ^+0^) — ["-(k)!; ^^^ the pairing gap in Eq. ([M)) 
is either Ami or Am3- This result does not rely on the 
approximation made in Eq. (|64p . because along this line 
Xk = or TT, and, as is clear from the structure of the 
matrix Eq. (|5n)) . the value of the gap along this line is 
indeed Ami or Am3- 

We plot the variation of f{6) ~ — 7i_(k)/77+(k) along 
the electron Fermi surface in Fig. [15] for the fitting pa- 
rameters given in Table IIXI Along the zone diagonal this 
expression indeed equals —1 in both cases; this means 
that the value of the gap is Am3 • 



it-jkp (9)) 

-ii+(kf (ft)) 

1.0 




-1.0 



u_lkp (ft)) 
-u+iki: (ft)) 



a) 




FIG. 15: The gap on the electron Fermi surface is given 
by Eq. (|64|l . Here we plot its angular dependent part, i.e., 
— 77_(k)/u+(k), at the Fermi surface, where i\^ =0. 



Fig. [TSJshow that f{0) is monotonically increasing as 
the angle 9 is varied from to 7r/2. The extremal values 
of the superconducting gap along the Fermi surface are 
therefore at 9 = and 9 = tt/2. As mentioned, /(O) = 
— 1, therefore one extremal value for the gap is Am3- The 



•^nodal 



AM3[meV] 
20[ 




AMi[meV] 



•^nodal 
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We compare these results to the ones obtained without 
the approximation leading to Eq. ([M)) . As shown in Fig. 
[TBI the transition between the nodal and the nodeless s- 
wave superconductor is indeed well captured by a straight 
line with a slope of —4.14 in agreement with the approx- 
imation in Eq. (|64p . The four superconducting phases in 
the plot do not break any space group symmetries be- 
cause they all correspond to Aig pairing. While at any 
finite temperature there is no phase boundary separating 
them, at T = the ground state energy is non-analytic 
at the transition lines, implying that these are distinct 
quantum phases. 



B. Superconductivity and spin-orbit coupling 



FIG. 16: The zero-temperature phase diagram for the super- 
conducting phase of iron-pnictides. Two parameters. Ami 
and Am3 determine the superconducting gap. The horizon- 
tal transition line is given by Am3 ~ 0; the vertical line is 
determined numerically, and for Ams J; lOmeV it coincides 
with the line where the zeroth order expression for the gap 
in Eq. (|63|l vanishes. In this figure we chose Ar > 0. When 
Ar < 0, only the sign of the gap on the hole Fermi surfaces 
changes, hence regions with s++ and s+_ gap exchange places 
in the phase diagram. 



other extremal gap value, at 9 = 7r/2, is 

A(V2) = -(Ami + Ams) + 2(^mi - AM3)/(7r/2). 

(65) 

From these two gap values we can determine the nature 
of the s-wave superconducting state. Let us assume, for 
concreteness, that the gap on the hole Fermi surfaces is 
positive, Ar > 0. If both Ams and A(7r/2) are posi- 
tive, then the gap on the entire electron Fermi surface is 
positive. This corresponds to the s_|_+-wave state. Anal- 
ogously, if both Am3 and A(7r/2) are negative, then the 
gap is negative on the entire electron Fermi surface, and 

this is the s^ -wave state. Lastly, when the gaps at their 

extremal values have different signs, the gap value, be- 
ing real, must vanish somewhere along the Fermi surface 
— the spectrum of such an s-wave superconductor has 
nodes. 

The transitions between the nodal and nodeless s-wave 
states occur when either Am3 or A(7r/2) vanish. The 
latter occurs if 



iM3 



/A 



Ml 



./(^/2) + 1 
./(V2)-f 



(66) 



For 6* 7^ and 9 ^ tt, the value of f(9) is model de- 
pendent. We find /(7r/2) = 0.611 (Ref. [ll|) or /(7r/2) = 
0.067 (Ref. \42t) . which implies the that the ratio in Eq. 
(|551) is —4.14 and —1.14, respectively. 



In the previous subsection we considered only spin sin- 
glet pairing terms with a focus on the s-wave supercon- 
ductivity. In this section we turn our attention to spin- 
orbit interaction which necessarily induces spin triplet 
pairing terms. In the absence of any spin-orbit inter- 
action, there is, of course, spin SU(2) symmetry, and 
therefore the spin triplet pairing terms are not induced 
by the spin singlet pairing. Once the spin-orbit interac- 
tion is included, the spin SU(2) symmetry is lost, and 
the spin part of any spin triplet pairing term transforms 
together with its orbital part under the operations of the 
space group. Therefore, some spin singlet and some spin 
triplet pairing terms may belong to a same irreducible 
representation of the space group. When this is the case, 
the presence of a spin singlet pairing, necessarily induces 
all spin triplet pairings which belong to the same irre- 
ducible representation as the spin singlet. This is concep- 
tually similar to the effect of the spin-orbit interaction on 
a spin-density wave order discussed in the section |lVl in 
the absence of the spin-orbit interaction, two spin-density 
wave orders with mutually perpendicular magnetizations 
are never related to each other, however, once the spin- 
orbit interaction is present, the two spin-density waves 
may belong to the same irreducible representation and 
necessarily induce each other. 

We begin by constructing spin triplet pairing bilinears 



H, 



triplet 



k 



<(-k)(icr2^)a0 



Mabi^'fsik) 



(67) 



where the three components of Mat are antisymmetric 
matrices, and fAaC^), with a = 1,...6 are the compo- 
nents of the "spinor" in Eq. (|17p with the spin label a. 
Unlike in the previous subsection, where the spin singlet 
was invariant under all space group operations, here we 
have to ensure that the spin part of the bilinears in Eq. 
(|67p also transforms under space group operations. The 
equation analogous to Eq. (|46l) for spin triplet pairing 
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terms is 



E^-,'-^(^) 



cdj 



EC(-k)(za2a,)„^V^(k) 



Y,D^nR)D^,!iR) 



cd 



E^;(-k)(za2a.)„^^,'(k) 



(68) 



The spin part of the triplet pairing term in Eq. (|68p trans- 
forms as an axial vector, i.e., 



E^-(^)<(^) 



E^^(-k)(*a2CT,)„^V'5^(k) 

k 

EfE^-(^)^S(-k)j La^Y^Df^ima: 
X (E^M(i?)V''(k) ) . 



Q/3 

(69) 



To show this, we use the transformation property of the 
spinor under a mirror reflection through the plane per- 
pendicular to a unit vector m: 

i?:f (a-)VMk) - ^a6(a-)(*m . aU4{k'), (70) 

where k' is the result of the mirror operation a"^ on k. 
Substituting this into the first line of Eq. (|M)) . we find 
that the spin part of the bilinear in Eq. (|68p transforms 
as 



{i(T2(Jj) — > (im • a) {ia2(Jj){im ■ (?) 



(71) 



under the mirror cr™. Because (cr;)-^o'2 = ~<72<^i for any 
I = 1,2, 3, we have 



{i(J2(Tj) — > i(J2{rrv ■ a)ai{vs\ ■ a). 



(72) 



Hence, the component of Mab which is parallel to the 
m is invariant under the mirror refiection cr™; the other 
two components of M°'^, which are parallel to the mirror 
plane, change sign. This is a property of an axial vector. 
Similarly, by replacing aj with 1 in Eq. (j68p . it is readily 
seen that the bilinears in Eq. (|46l) are spin singlets. 

The second line in Eq. (|69|) tells us that the symme- 
try properties of a bilinear are given by the product of 
the irreducible representations of the orbital part of the 
bilinear and the axial vector representation of the spin 
triplet. The symmetry properties of the orbital part of 
pairing bilinears have been presented in the previous sub- 
section. The symmetry properties of the spin part, an 
axial vector, are given by the A2g, when the spin points 
in the z-direction, and by the Eg, when the spin is in X- 
or y-direction. 



1. Superconductivity and spin-orbit coupling nearV 

At the F-point only the orbital A2g bilinear is antisym- 
metric in the (orbital) components. When Mat points in 
the z-direction, the entire spin triplet pairing term at F 
transforms according to A2g (8) ^2g = ^ig- When Mat 
is perpendicular to the z-axis, the spin triplet pairing 
term transforms according to Eg ^ A2g = Eg. We are 
interested in the former, i.e., the Aig, since such a pair- 
ing term has the same overall symmetry properties as 
the spin singlet s-wave considered earlier. The resulting 
pairing Hamiltonian reads 



ffr,tripict = ^rt 2^ V'r4(-k 



-i 
i 



j ipr.^Ck.) + h.c. 
(73) 



where by time reversal symmetry, which we assume to be 
present, Art is real. 

We can use the same Nambu spinor for the states near 
the F-point as in Eq. ()54p . The Bogolyubov-de Gennes 
Hamiltonian is 



H, 



r,BdG — 



/ir(k) + iArT2 Arl + ArtT2 



(74) 
where T2 acts in the same space as /ir(k) defined in Eq. 



1.0 




FIG. 17: The angular dependence of Ar/2ylk along the inner 
(top) and the outer (bottom) hole Fermi surfaces, calculated 
for the parameters from Ref. [ill . Table IIXI The variation of 
this ratio enters into the formula for the anisotropy of the 
superconducting gap on the hole Fermi surfaces, Eq. (|79p. 

Using the shorthand notation 

/ir(k) + iArT2 = Akl + f-Sk, 



Ar]L + ArtT2 = Ck + f-Dk, 



(75) 



we use the Eq. (jAip in the appendix. The eigenvalues 
of the Bogolyubov-deGennes Hamiltonian Eq. ((7^ are 
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given by Eq. (|M1) . 



El =Al + Bl + Cl + Dl 



±2 



(v4k-Bk 



Ck-Dk 



BlDl- 



[b^-d^' 



(76) 

The normal state dispersions for the two hole bands in 
the presence of the spin-orbit coupling are 



.(±) 



Ak ± \Bv 



Av±Bv 



(77) 



where we write i?k = |-Bk|. We can therefore elimi- 
nate Bk in the Eq. ((76)) in favor of ±(e{^ — A\^), while 
noting that the scalar product B\^ ■ D^^ — ^Xr^rt is k- 

independent. At the Fermi surface, one of the ej^. s van- 
ishes, and the gap is given by 



aFS 



2Al 



Al 






K 



^IMl 



K 



iA2)-|-^kArArAri 



1/2 



(78) 

In the above expression, k lies either on the smaller or 
on the larger hole Fermi surface, which are not circular. 
Because Ak = er + k^/(2mr), is circularly symmetric, 
the gap on cither of the Fermi surfaces is anisotropic. In 
the limit when Apt vanishes, however, the gap is isotropic 
and given by Ar , even for a non-zero spin-orbit coupling 
constant Ar- 

To determine the size of the anisotropy, we note that, 
for the parameters of Ref. [ill, and for Ar = 80meV, the 
magnitude of Ak is at least 40meV on both hole Fermi 
surfaces. Therefore, for A's small compared to this en- 
ergy scale, we can expand Eq. ([78l) in powers of Ar/Ak, 
and Art/Ak. In this limit we find that the gap on the 
Fermi surfaces is given by 



aFS 



At 



2Au 



-A 



rt 



(79) 



In Fig. (|T7)) we plot Ar/2Ak on the inner and the outer 
hole Fermi surfaces for the parameters of the model in 
Ref. [ll| given in Table IIXI This approximate expression 
vanishes if Ar/Art intersects either the upper blue curve 
or the lower red curve in Fig.[T71 Therefore, to this order 
we would conclude that either the inner or the outer hole 
Fermi surface has gap nodes. Closer analysis of the ex- 
pression Eq. ([75]) reveals that the true gap nodes appear 
only under much more stringent conditions: Ar = Art, 
and Ak = Ar/2 at some k on a Fermi surface. There- 
fore, the gap nodes determined from the approximate 
expression Eq. (|79p are only "near gap nodes", in that 
". . ." terms give the gap a finite but small value (of or- 
der Ap/A^). Nevertheless, because Ar/^k has opposite 
sign on the inner and the outer hole Fermi surface, we 
see that, if the "near gap nodes" appear, then they only 
appear on one of the Fermi surfaces; the gap on the other 
Fermi surface remains relatively isotropic. 



2. Superconductivity and spin-orbit coupling near M 

At the M-point we have to consider the direct prod- 
ucts: £'mi (8) i?Mi, -Ems <^ -Ems, and i?Mi (8) -Ems- These 
are given in Table IIVI We see that there the only prod- 
uct which contains Eg is E'mi <8i i^Ms- Hence we can 
combine it with the planar. Eg, components of the ax- 
ial vector Mab and form an Aig pairing bilinear. Since 
it is also antisymmetric in its orbital indices, it satisfies 
the proper anticommutation relations, and such term will 
necessarily get induced by the spin-orbit coupling in the 
presence of the Aig singlet. Therefore, the spin triplet 
pairing Hamiltonian at the M-point reads 



-f^M, triplet =A]V[t 2^ [ ^X4.(^l^) 



-i 
-1 



V'y;(k) 



■V'yt(-k) 



1 

-i 



i'xtik) + h.c 



(80) 



The triplet pairing parameter Am* is real by time rever- 
sal symmetry, which we assume to be present. 

We now construct the eight component Nambu spinor 



*M(k) 



/ ^xt(k) \ 

V'y;(k) 

V^xi(-k) 



(81) 



which we use to set up the Bogolyubov-dcGcnnes Hamil- 
tonian 



ffBdG,M = X!*M(k) 

k 

where 



n{k) 



-n{k) 



*M(k), (82) 



Ami -iAm* 

Am3 -Am* 

-Am* Ami 

iAm* Am3 



n{k) 

A 



h+,{k) A 
At h^{k) 



2^^ I 1 



(83) 

(84) 
(85) 



Because the matrix in Eq. (|82|) can be writen as t^ (g) 
■H(k) + Ti (g) A. it manifestly anticommutes with r2 (g) 1. 
Therefore, for any eigenstate at k with energy E^, there 
is an eigenstate at k with energy —E^. 

To proceed with the approximate solution to the spec- 
trum, we first perform a unitary transformation which 
diagonalizes h^Ck) and h^{k); for the parameters of 
our model, the two eigenvalues for each are split by an 
energy which is larger than both the spin-orbit coupling 
energy scale and the pairing scale. Therefore, we can 
safely project onto the states which are near the Fermi 
level. This will result in a 4 x 4 Hamiltonian with the 
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particle hole symmetry, meaning that the secular equa- 
tion will reduce to quadratic. The corresponding basis 





(86) 



into Eq. (|M|I . 



El=Al 



Bi + a 



i + Dl 



±2^/(AkBk 



C^D^ 



BlDl 



BvDy 



(97) 



where /i]y[(k)|±) = e^. |±). In this basis, the effective 
Bogolyubov-de Gennes Hamiltonian reads 



H 



where 



M3dG 



H(k) = 



H(k) A(k) \ 

A(k) -H(k); 






(87) 



2AMK C|j. 



~(1-) 



-AMtA^* A(i-) 



(89) 



where A*^^"*") was given in Eq. (|64p . and analogously, 

A(i~) =^(Ami + Ams) + ;j(Ami - Am3)-t]- ^" 



r(i-) 



u,+ (k) 
(90) 



where ui's are defined as h (k) = z/;+(k)l + w_(k)(7^ 
w_(k)cr2'; 



i 

1 



(91) 



The phase of this complex number may be eliminated 
by a unitary transformation generated by \® gj,. This 
transformation leaves invariant all terms which do not 
contain k\ in all other terms k is transformed into 



.(k) 



w_(k) 



1/2 



(92) 



6f+)-zi+(k)ef-)-u;+(k), 



The Bogolyubov-de Gennes Hamiltonian in Eq. 
can be written in the form tj, (g) ■H(k) +t\® A(k). Eigen- 
values of such matrices are derived in the Appendix [A] 
The spectrum of Eq. ([5^ can be obtained exactly by 
substituting 



:i |';(1 + ) 



A=-(ek "+f'k 



(1-) 



Sk=^(AM|«|,0,ef+) 



r(i-) 



i3k = ('-AMt|^c|,0,i(A(i+)-A(i- 



(93) 
(94) 
(95) 
(96) 



~(i+) 

^k - "^k 



The dispersion in the normal state is e^ = Ak± |i?k| — 
Ak ± -Bk, where we, just like in the previous subsection, 
define i?k = |-Bk|. Because along the zone diagonal 
ej^ , it is the non-zero spin-orbit coupling Am 
that lifts the crossings between the two "elliptical" elec- 
tron Fermi surfaces. The "-I-" sign therefore corresponds 
to the inner electron Fermi surface, and the "— " sign 
to the outer Fermi surface. Eliminating A\^ in favor of 
e^ ^ _Bk and noting that at an electron Fermi surface, 
one of e^k vanishes, we find that the gap is given by 



aFS _ 



2Bl + Cl 



Dl- 



2\/5k + Bl[Dl- (bk • D^Y T 2Ckbk • D^ 



ClDl 



1/2 



(98) 



where we defined the unit vector bk = i?k/-B + k. In this 
expression, the momentum k lies on one of the electron 
Fermi surfaces. 

The magnitude of Sk is at least 28meV on cither of 
the Fermi surfaces for the parameters in Table IIXI cor- 
responding to Ref. [m and Am = 80meV. For A's small 
compared to this energy scale we can expand the Eq. ([SS]) 
in powers of AMi/^k, AMsZ-Bk, and AMt/-Bk- In this 
limit we find that the gap on Fermi surfaces is 



aFS 



= |Ck±bk-i?k| 

= |aMi(k)AMi - 



aM3(k)AM3 + aMt(k)AMt| 



(99) 



Here ". . ." represents terms of order A^/Bk- 

The second line of Eq. ([99|) defines the coefficients 
aj (k) , which are plotted in Fig. [18] as a function of an- 
gle 6 on the two electron Fermi surfaces: inner (blue) 
and outer (red). For comparison, in Figs. [T8b andlTSb. 
we also plot Uj 's in the absence of spin-orbit interaction 
(dashed curves). The most important consequences of 
the spin-orbit coupling on the superconducting quasipar- 
ticlc spectrum are: a) coupling of the two electron Fermi 
surfaces via lifting of the degeneracy which occurred for 
Am = at the Brillouin zone edge — the Fermi surfaces 
reconstruct and each Fermi surface is four-fold symmet- 
ric; b) the admixture of the spin triplet pairing. Am*, 
which contributes to the overall size of gap with a co- 
efficient QfMt- and which is comparable to the two spin 
singlet coefficients q;mi, Q;m3 in the Eq. 
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c) 



FIG. 18: The coefficients a(k) in the electron Fermi surfaces 
gap estimate, Eq. (|99p . as a function of angle 6. The coeffi- 
cients imply a four-fold symmetric gap on each electron Fermi 
surface. The red and blue line correspond to the outer and 
inner Fermi surfaces respectively which do not touch in the 
presence of spin-orbit coupling. The dashed lines in a) and b) 
correspond to no spin-orbit interaction, in which case OfMt is 
absent. 



VI. DISCUSSION 

We applied the method of invariant a^^i^^ to construct 
an effective low-energy effective continuum model for iron 
based superconductors. Instead of repeating the conclu- 



sions which were stated in the introduction, in this sec- 
tion we discuss directions along which this line of rea- 
soning can be extended. As we have mentioned, among 
the advantages of our approach is that it readily allows 
one to analyze the effects of an external magnetic field. 
Therefore, unlike previous continuum models, ours has 
the advantage of being suitable for a calculation of cy- 
clotron resonances. 

While for the spin-density wave order we have mapped 
out all of the symmetry distinguishable order parameters 
with wave-vector M, we focused primarily on the Aig 
superconducting order which preserves the time-reversal 
symmetry. In the Appendix |D] we list the remaining 
symmetry distinguishable pairing terms without analyz- 
ing their effect on the quasi-particle spectrum. Since in 
some of the iron based materials there is a possibility of 
having a non-Aig superconducting order, such analysis 
can be readily performed using our low-energy effective 
theory. 

The method of invariants used here is generioi^. It re- 
lies on identifying the symmetry of the exact Bloch states 
at isolated points in the Brillouin zone, and an expan- 
sion in powers of momentum away from such points. As 
long as the Fermi momenta arc a small fraction of the 
Brillouin zone, low order expansions are sufficient. We 
demonstrated the applicability of the method by fitting 
to two ti ght -binding band structure calculations reported 
in Refs.llll (8 bands) and IH (5 bands). The method 
can be applied to construct a low-energy effective theory 
based on a microscopic band structure for any other iron 
based material with the same lattice symmetry, but pos- 
sibly different spectrum with small pockets. To do so, 
one simply has to select the isolated points in the Bril- 
louin zone and the corresponding states near the Fermi 
level as the basis states. Then project the Hamiltonian 
onto this basis, and perform the expansion in powers of 
momentum away from such points. If some states do not 
cross the Fermi level, but lie close to it, one could include 
those as well by extending the number of components of 
our "spinor" ■0k, cr, Eq. ([T7|) . 

The space group wc study, PA/nmm, is the group of 
symmetries for the parent compounds with a single iron 
layer per unit cell, such as those within the 1111, 111, 
and 11 families^. Since this group is non-symmorphic, 
the conclusions we draw for these compounds apply di- 
rectly. For the 122 family the space group of symmetries 
is lA/mmm. The puckered iron-pnictide layer is still the 
basic structural unit, however, the puckering alternates 
along the z-direction from layer to layer. As a result, 
this group is symmorphic, and therefore our results must 
be modified when applied to such materials. In the limit 
when the inter-plane coupling can be ignored, our results 
apply without modifications. While we have implicitly 
assumed absence of /cz-dispersion, our low-energy effec- 
tive theory has the same form, with possibly modified 
parameters at fc^ = and k^ = tt/c. For weak inter-layer 
coupling, as applicable to these systems, the k^ = and 
kz = """/c planes in the Brillouin zone should be folded. 
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Our results can then be used as an appropriate starting 
point, defining the basis which is weakly mixed. 

In the paper we present the irreducible representations 
of PA/nmm the way we obtained them. We followed the 
physically motivated method of C. Herring^, by con- 
structing the factor groups for non-symmorphic space 
groups, and the method of induced representations out- 
lined in Ref. Q. Alternatively, using different methods, 
generators, product tables and irreducible representation 
for all 230 space groups — including PA/nmm — along 
the high symmetr y li nes, or at high symmetry points, 
appear in the Ref. [Sll Our results are fully consistent. 

The Landau free energy function, built based on the 
irreducible representations at the M-point, displays two 
distinct isolated extrema below the ordering temperature 
Tm. One corresponds to the collinear spin-density wave, 
in our notation transforming as i^^i ^r ^-mij whose 
presence in iron-pnictides have been well established^i^. 
The other corresponds to the coplanar, four-fold symmet- 
ric, spin-density wave order, E^^±E^^, which may have 
been observed recently^. When the former order pa- 
rameter corresponds to a minimum, the latter is a saddle 
point of the Landau free energy, and vice versa. Which is 
preferred depends on the sign of the coefficient A' in the 

quartic term. A' I A^x ~ A^y ) . We speculate that 

the transition observed in Ref. |35| is due to the tempera- 
ture dependence of A' which is negative above the copla- 
nar ordering temperature and positive below. Such tran- 
sition is necessarily first order, in agreement with Refs. 
ISSl and |40|, although without invoking any spin-nematic 
model. Our method allows us to study the difference in 
the electronic spectrum for the collinear and the coplanar 
spin-density waves. Interestingly, while the former breaks 
the time-reversal symmetry, the electronic spectrum is 
still Kramers degenerate. That is because, the product 
of time-reversal and inversion is preserved. On the other 
hand, the coplanar spin-density wave breaks both the 
time reversal and inversion, as well as their product. The 
resulting electronic spectrum has no Kramers degenera- 
cies. Therefore, in this phase all Fermi surfaces split as 
seen in Fig. [TH We expect this to have interesting ram- 
ifications for the nature of the superconducting state if 
it microscopically coexists with the coplanar spin-density 
wave. 

Another advantage of the low-energy effective model 
is that it is able to capture the anisotropy of the s- 
wave superconducting state with a minimal number of 
]i.-mdependent parameters. The anisotropy of the quasi- 
particle gap on the electron pockets, and consequently 
their nodal structure, depends on the ratio of Ami and 
Ams, as seen in Fig. 1161 This holds even when the spin- 
orbit coupling is ignored. The spin-orbit coupling intro- 
duces an additional parameter. Am*, corresponding to 
the pairing strength in the spin triplet channel which is 
necessarily induced. The spin-orbit also removes the de- 
generacies at the Fermi surface crossings, thereby split- 
ting them. The resulting quasi-particle gap is therefore 



unambiguously four-fold periodic, see Fig. [181 A similar 
conclusion was obtained in Ref. |5J where a phenomeno- 
logical spin-independent hybridization parameter, X{(p), 
was included. We see that, in the absence of spin-orbit 
coupling, the two Fermi surfaces are hybridized only if 
the PA/nmm lattice symmetry is broken. In the presence 
of such symmetry, the hybridization parameter necessar- 
ily arises from the spin-orbit coupling, and its detailed 
spin-structure is captured within our approach. 

On the hole pockets, a single Aig spin singlet pairing 
parameter, Ar, yields isotropic quasi-particle gap. In the 
presence of the spin-orbit coupling, the spin triplet, k- 
independent, pairing, parameterized by Art, is induced. 
The addition of this term leads to the gap anisotropy 
on the hole Fermi surfaces. In this case, we find the 
anisotropy to have a different strength on the two hole 
pockets. This result may be tested using spin polarized 
ARPES. 
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Appendix A: Eigenvalues for Bogolyubov-deGennes 
Hamiltonians 



In Section |Vl the spectrum of superconducting states 
was determined from the eigenvalues of the Bogolyubov- 
deGennes Hamiltonian. Once spin-orbit coupling is in- 
cluded at either the F point or the M point, such Hamil- 
tonian matrix can be casted in the form 



i?BdG = 5I*^W 



T-3 ® Uk + Sk • CT 



+ri® (Ck + ^k-<?)]*(k), (Al) 



where for each k, Ak and Ck are scalars, while B^ 
and Dk are vectors "dotted" into the three Pauli ma- 
trices a = (ci, (T2, 173). Here we demonstrate a simple 
method for obtaining the eigenvalues for such a 4x4 ma- 
trix Hamiltonian. 

Because the Eq. (jAip anticommutes with r2 1, for 
each eigenstate at k with an eigenvalue i?k, there is an- 
other eigenstate at k with the eigenvalue —E\^. There- 
fore, solving the secular polynomial reduces to finding 
zeros of a quadratic function. 

Squaring the matrix in Eq. (|A1[) . and using the well- 
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known commutation relation [(Ta,(T(,] = 2ieabc'^c, gives 



H^ ^U[Al + Bl + Cl + Di 



t®a-[ 2AkSk + 2Ck£'k \ 
-2t2®S- (^Bk X I?k j ■ (A2) 

Choosing a basis in which T2 is diagonal, the eigenvalues 
of the matrix in Eq. (|A2[) . £'^, can be readily read off. 
Because i3k x ^k is orthogonal to both B\^ and to I?k, 
the result is independent of whether the eigenvalue of T2 
is +1 or —1. We thus find 

± 2^ (A^k + Ck^k) ' + BlDl - (Sk • ^k) '. 

(A3) 

Therefore, the four eigenvalues are ± \J E^ with uncorre- 
lated ± signs inside, and in front of, the overall square 
root. 



Appendix B: Group Pm and the construction of its 
irreducible representations 

Since the space group of iron-pnictides lattice symme- 
tries, P4/nmTO, is non-symmorphic, there is a particular 
procedure, due to Herring, for the construction of the 
space group irreducible representations at the Brillouin 
zone boundarji^. Wc follow the procedure as illustrated 
in Ref. [M 

The main idea is to take the subgroup of 'even' trans- 
lations, 

rM = {{e|t}|exp(zM.t) = l}, (Bl) 

and decompose the space group into cosets 

Q = RiTm + R2TM + ■■■■ (B2) 

Since 7m is an invariant group of Q, the cosets form fac- 
tor group Pm = G/Tm- The odd translation, {ejlO}, 
is among the elements of Pm, and for any irreducible 
representation at M, this symmetry operation must be 
represented by 



D({e|10}) 



iaM.-x 



-t. 



(B3) 



The irreducible representations at the M-point are there- 
fore all irreducible representations of Pm for which Eq. 
(jB3p holds. All the other irreducible representations of 
Pm are unphysical and should be disregarded. Our task 
here is to identify group Pm and select its physical irre- 
ducible representations. 

The group Pm has 32 elements; in Seitz notation 
these are 16 operations, {^Jt}, already present in Pp, 



and an additional 16 operations obtained by multiply- 
ing the 16 elements of Pr by an 'odd' translation: 
{e|10}{g|T} = {g\10 + t}. Note that the multipli- 
cations rules are not the same as in Pr, as the two 
symmetry operations {giti} and {5|t2} are equivalent 
in Pm only if exp (iM • (ti — 12)) = 1. For example, 
{cr^|00}{i|ii} = {Cf |ii}. Wc found that Pm is not 
isomorphic to any crystallographic point group in three 
dimensions^. Therefore we were left with the task of 
deriving the irreducible representations for Pm- 

In constructing irreducible representations of an un- 
known group, a good head start is gained by finding the 
largest possible invariant subgroup whose irreducible rep- 
resentations we already know. In the case of Pm , wc no- 
tice that four elements of the group, {c^^l^^}, {f^lj^}, 
{cr^|ii}, and {e|10} commute with each other, and at 
the same time, each of these elements combined with the 
identity operation forms a cyclic subgroup isomorphic to 
C2. Therefore, a set of operations obtained by succes- 
sive action of these four elements forms a group which is 
a direct product of four cyclic groups, H = (C2)®'*, and 
is an invariant subgroup of Q. This subgroup has 16 ele- 
ments and is Abelian. Each cyclic group C2 has two one- 
dimensional irreducible representations, one even and one 
odd. The irreducible representations of Ji are then ob- 
tained by multiplying the irreducible representations of 
each cyclic subgroup. There are 16 one-dimensional rep- 
resentations given by 



^^«'"^^HK^iHrK'iH}"KiHne|ior) 



: f "7/"C^T^ 



(B4) 



where each ^, 77, ^, and r can take either -1-1 or —1 value 
corresponding to a representation that is even/odd under 

{c^^lH}' i'^^lll}-- i'^^lll}-- and {e|10}, respectively. 
Here, m, n, p, and q are integers, although it is sufficient 
to consider values and 1 to find all the elements of H. 
For brevity, we use only the sign, -I- or — , and drop I's in 
the designation of the irreducible representations of H . 

The next step is finding the induced representations. 
Wc write Q ~ RiH + R2'H, where Ri — {e|00}, while we 
choose i?2 = {cr^|00} for simplicity. Taking an irreducible 
representation of H, D^?'''^'^), the induced representation, 
B(ivCr) I g ig defined as the 2x2 matrix 



(^D^ini-^tg)^J{p) 



0, 



R7 PR, 



Ri'PR, i n ' 
(B5) 



where P £ Q. The induced representations are, in gen- 
eral, reducible. 

To determine whether given £)(^''?'^) t ^ is reducible or 
not, for each R^, we create an irreducible representation, 
i:)(?'?CT),(0^ defined as 



£,(«'7Cr).(»)(5) ^ £)(C'Kr) (^R-^SR^) , 



(B6) 
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with S G 7i. These irreducible representations are 
then separated into classes of equivalent representations. 
However, since in our problem there are only two irre- 
ducible representations created this way, we can either 
have the two of them being equivalent and belonging to 
the same class, or being inequivalent and each defining 
its own class. If the two irreducible representations de- 
fined in Eq. (JB6| arc inequivalent, Df?''?"^' 'f tj is an irre- 
ducible 2-dimensional representation of Q and we name 
it ^(«''C^). Otherwise, I?(«''C^) f 5 is reducible. In such 
a case, the irreducible representations are constructed 
by finding the unitary matrix U which connects the two 
equivalent irreducible representations, Ai and A2, 



£,(«rKT),(2)(5.) ^ [/-lL»(«'";^)'(l)(5)t/. 



(B7) 



Since all £'^^''^'^)'s are one-dimensional, the unitary ma- 
trix is just a scalar of modulus unity U = exp(i(^), and 
we choose it to be 1. The two irreducible representations 
obtained from £)(?''?'^) arc then given by 



For our problem we need first to find how {cr^jOO} con- 
jugates the elements of "H. Using the Scitz product rule 
in Pm, we find 



K|00}-Ha^|ii}K|00} =W'^\l^,} = {e|10}{a^|ii}, 

(BIO) 

K|00}-i{a^|ii}K|00} ={a^lii} - {e|10}{a^'|ii}, 

(Bll) 

KIOOI-^KIHIKIOO} =K|Ii} = {e|10}{a^|ii}, 

(B12) 



{f7"|00}-He|10}{fT"|00} ={e|10}. 



(B13) 



a[^'''^'\s) = 


D^^"'^^\S), 




a[^'''^^\r,s) = 


+i?(«"«"'(5). 


(B8) 


4«"C-)(5) = 


^(?r,Cr)(5)^ 




Ai^^'^^^R^S) - 


~D^^'"^^\S). 


(B9) 

1 



To derive these it is useful to notice that {cr^|00} ^ — 
{cr^ |00}. For an arbitrary element of H we find 






D^^^'i^^ (({ellOlKlii})" ({e|10}{a^|ii})" ({e|10}K|ii})'' {e|10r 



,y| 1 1 
I2 2 



_ fn m^p^m+n+p+q 



(B14) 



This equation defines irreducible representation 
£)(5';C'r),(2)_ This representation is equivalent to 

D{ivCr),ii) = oiinCr)^ provided that ^■^r^^(T>T't = 
^n^rrii^p^m+n+p+q £qj. ^^^ .^^ ^^ ^^ ^^^ ^ integer. Arrang- 
ing this equation and bearing in mind that ^, 77, C, and 
T only take values ±1, we find that the two irreducible 
representations are equivalent if 



i^Vr) 



m-]-n 'o 

r 



const, 



(B15) 



for any integer m, n, and p (also q, but it drops out 
of the condition) . The left hand side of Eq. (JBISP is p- 
indcpendent only if r = -1-1. Similarly, it is m and n 
independent provided that ^ry = +1. There are, there- 
fore, four irreducible representations of "H that each in- 
duce two one-dimensional irreducible representations of 



D(++-+)^g) = 4++-+) ©4++-+), 



D^-++^tg) = 4-~++)®4-++\ 



D( 



-+) 



to 



a[—+'^ 3) Ai,—+\ (B16) 



These irreducible representations of Pm sue defined as 



l^iivCr) 



({^""iHrKiHrKi^n^iior) 






A 



iivCr) 



= f "ry",^Pr«, (B17) 
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where the only allowed combinations of ^, 77, C, and t are 
those present in Eqs. (|B16p . 

The remaining twelve irreducible representations of % 
induce two-dimensional irreducible representations of Q. 
These come in six pairs of equivalent irreducible repre- 
sentations: 

Ei++—) ^ E^—+-\ S(+++-) ^ S( ), 

Ei+-++) ^ £(-+++), S(+— +) - £(-+-+)(B18) 

The matrices for the induced representations, i.e., the 
representations in Eq. (|B18[) . are 

ewx^r^ \ 

i?(«''?-)(Kioo}{cT^|iir{a^|iirK|iir{e|ior) 

Q pn m^p^m+n+p+q 

^m^uQ^^q Q 



(B20) 



The characters are 

X 



71 in ^p rn~\-n-\'p-\-q 



^"i^n^-p^q ^ ^n^m^p^ 



^m^n^P^g 



1 + (e^r)™+" tP 



(B21) 



A^vC-r) 



({a-|00}{a^|iir{cT^|iirK|iir{e|10r) = 0. 

(B22) 



Hereby we have finished the construction of the IR's for 
the unknown group Pm- 

From Eqs. (jB17|) . it follows that for all one-dimensional 
irreducible representations. 



4«''^^)({e|10}) = 4«'"^^)({e|10}) 



-1. 



(B23) 



Hence, all one-dimensional irreducible representations of 
Pm are unphysical. Similarly, we use Eqs. (|B19|) - (|B20[) 
to determine the representation for the odd translation 
to be 



^(?'(Cr)(|g|iO}) = (-l)^l 



(B24) 



in two-dimensional irreducible representations in Eqs. 
(jBlSp . Therefore, the irreducible representations for 
which T = —1, i.e., those in the last line of Eq. (jBlSp . 
are unphysical. The remaining four irreducible represen- 
tations in Eqs. (|B18p arc physical and we relabel them 



Em3 = E^++- 



Em2 = E'^^ 

^M4 = E^ 



(B25) 



Let us just note here that the first member of each E'm 
doublet is odd, while the second member is always even 
under {<J^\h\}- There are four combinations for the par- 
ity under {cr^|^^} and {cr^|ii} and they are all present 



in Eq. (jB25p . Per construction, all the physical irre- 
ducible representations at the M-point are odd under 
{e|10}. Finally, the {cr^ |00} mirror is represented by ci 
Pauli matrix in each i^M irreducible representations, i.e., 
it swaps the two components. 

The irreducible representation table, Table Hill as well 
as the irreducible representation multiplication tables. 
Tables HVl and FVl in the main text, follow directly from 
Eqs. (|BT9)) and ((B20| . 

The unphysical irreducible representations we con- 
struct here are, in fact, not entirely unrelated to the space 
group. Noticing that for all of them, £)({e|10}) = +1, 
it follows that these irreducible representations are phys- 
ical at the F-point. Each irreducible representation of 
D4/1 is, therefore, equivalent to one of the (unphysical at 
M) irreducible representations of Pm derived here: 



(+H 






A\ 
A 



B 
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B. 



2s 



E,=E 



(++++ 
2 

{+— + 



A. - A^ +' 

B2u = A[ 

E,. = E^-+++ 



+) 



(B26) 



Appendix C: Particle- hole bilinears in the 
low-energy effective theory 



i i.R.™s r^x^^ 



4 + 

Btu 

Et 
E+ 

En 

^M2 

P+ 

E+ 



(m) -p(m) 



Arl2, Am12, Amo"3 
Ar<7'2 

^1 v/2 

Arfi 

ArcTa, A3I2, A3(T3 
- 12-0-3 

^'^^ 

-Aig2 + A2<Ti Aio-2 + A2q' 
V2 ' V2 



-A3 + A]vi 



0-2 



-^3 + -''M 
V2 



0-2! 



( 



\ V2 ''^^ V2 "V 

[^*~7^' ^ 1 

\ l2~'^3 — A(3CT2 + A7q-i 

A4CT1+A5CT2 \ . I2-CT3 

V2 '^"5 V2 

-A4a2 + -^5<^1 \ 3i-2~fT3 

^^ '^^~72~^ 

' \ 3L2~l~g"3 AgiTi — A7g2 

/**^2~' Tf 

' \ l2 + °'3 AsCT2 + -''7gl 

'^'^"^2^' 75 — 

A4CT1— A5a-2 \^l2+£3 



V2 



V2 



A40-2 + -^501 \ ^2+0-3 



Here we list the matrices T^^ according to their sym- 
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mctry properties. The eight Aj's are "Gell-Mann" matri- 
ces defined as 



(CI) 
(C2) 
(C3) 
(C4) 




A4 = 



A5 == 



^0 
Xj= \0 -i 
.0 i 




1 
As = ^ I 1 
-2, 



We also used their combinations 



^M 



V3 



A8 + 



diag(l,l,0), 



(C5) 



Ar = -y ^As + -i=l3 = diag(0, 0, V2). (C6) 



Appendix D: The classification of all pairing terms 
in the low-energy effective theory 

In Section |Vl general forms for pairing bilinears in the 
low-energy effective theory are presented in Eqs. (|45p and 
(|67)) . The focus of that section is, however, on the Aig 
pairing terms. Here, we tabulate, according to their sym- 
metry properties, all other spin singlet pairing terms that 
may appear in the low-energy effective theory. 

A general spin singlet pairing term in our low-energy 
effective model has form 



^Vr(-k)A^V't(k), 



(Dl) 



The symmetry properties of the pairing term are deter- 
mined by the matrix 



M 



Mr 

-'^'PDW 



A^PDW 

Mm 



(D2) 



The 2x2 matrix Mr corresponds to the pairs formed 
by the hole band states. Its elements are classified ac- 
cording to 



E„ 



'En = A 



19' 



i2s 



B 



19 



.B. 



2g, 



(D3) 



where the orbital antisymmetric combination A2g should 
be dropped. Similarly, pairing of states in the electron 
bands is specified through the 4x4 matrix A^m- Its 
elements are classified according to 



Ml ^ Emi 


= Aig ® B2g (B A2u (B Biu, 


(D4) 


M3 <8l Em3 


= Aig ® B2g (S) Alu (S) B2u, 


(D5) 


Ml ^ Em3 


^ Eg®Eu, 


(D6) 



where we drop the antisymmetric orbital combinations, 
A2u, and Aiu in the first two equations respectively. 

The 2x4 matrix A^pdw corresponds to pairing terms 
between a hole and an electron state. Such pairs have 
momentum M, and their symmetry properties are given 
by one of the E-m. irreducible representations found in the 
decomposition 



Eg (8) i^Mi = Ems ® £'m4, 
En (8) Em3 = Emi ^ Em2- 



(D7) 

(D8) 



The pairing terms with no overall momentum, classi- 
fied according to their symmetry properties, are defined 
by the following M matrices: 



119 



B 



19 



"1 0' 
1 


' 


"10 0- 

10 

.0 0. 


' 


"0000 
10 


.0001 


"0 r 
1 


7 









s. 



29 



1 



Lo 



E„ 



A. 



B' 



2u 



E,, 





10 

10 



10 



10 





1 


10 





1 

10 







-1 







. 















1 







-1 







1 









(D9) 



(DIO) 





-1 

(Dll) 
(D12) 



(D13) 



(D14) 



(D15) 



Any 2x2 matrix should be substituted in Eq. (|D2p in 
place oiMr with the other block matrices being identical 
to zero; any 4x4 matrix in this list should be substituted 
in Eq. (|D2p in place of A^m with the other block matrices 
substituted by blocks of zeros. 

In Section |Vl wc discussed the Aig pairing terms, Eq. 
(|45)) . and their consequences on the low-energy effec- 
tive model quasi-particle spectrum. Similarly, one could 
construct i?2g pairing terms which would be described 
by three parameters, each multiplying one of the pair- 
ing terms following from Eq. (jDlip . The i?2g super- 
conducting order leads to a gap formation on both hole 
and electron Fermi surfaces of the low-energy effective 
model. The symmetry properties of B2g dictate that 
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quasi-particlc gap is odd under mirror reflection oper- 
ations {cr^jOO} and {cr^jOO}, implying that the gap on 
the hole Fermi surfaces has nodes for k^ ~ or ky = 0. 

For any other irreducible representation in the list of 
the spin singlet pairing terms with no overall momentum, 
there is always just a single term present. The supercon- 
ducting order corresponding to one of these irreducible 
representations would therefore produce a quasi-particle 
gap on either the hole or the electron bands, but it would 
never affect the spectrum of both hole and electron ex- 
citations. This suggests that pairing into any of these 
channels is likely irrelevant for the iron based supercon- 
ductors. 

Here we list the complete set of the symmetry adapted 
matrices jMpdw: 



Emi 

Em2 
Em3 
Em.4 



r 







' 


0] 




1 


' 



0' 




10 


' 


10] 







' 






0] 





-10 





-1 ol 









-1 








ol 























0' 


-1 












(D16) 
(D17) 
(D18) 
(D19) 



There is precisely one doublet of pair-density wave bilin- 
ears for each irreducible representation at the M-point. 
Therefore, any spin singlet pair-density wave supercon- 
ductivity in iron-based superconductors built upon the 
low-energy effective theory can be described by a single 
two-component k-independent order parameter. An on- 
site pair-density wave located on iron sites corresponds 
to one of the bilinears in Ec^. (|D19[) : similarly, an on- 
site pair-density wave on pnictide sites is represented by 
bilinears in Eq. (|D17p . 
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